TABELA 6 — Viga continua de dois tramos, parcialmente encastrada nos apoios extremos, com cargas 
uniformemente distribuidas 


la/ 1; 


1,0 


1,1 


Voa 


1,3 


1,4 


1,8 


1, 


20 


2,25 


2,5 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,1250 


0,1388 


0,1550 


0,1738 


0,1950 


0,2188 


0,2450 


0,2738 


0,3050 


0,3380 


0,3750 


0,4766 


0,5938 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0625 


0,0595 


0,0568 


0,0543 


0,0521 


0,0500 


0,0481 


0,0463 


0,0446 


0,0431 


0,0417 


0,0385 


0,0357 


0,0000 


0,0000 


0,0009 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0625 


0,0792 


0,0982 


0,1194 


0,1429 


0,1688 


0,1969 


0,2275 


0,2604 


0,2956 


0,0333 


0,4381 


0,5580 
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0,0000 


0,0009 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


0,0000 


TABELA 6 — Viga continua de dois tramos (cont.) 
n= 1/4 
la/h ES CIERETEaNaKAE: 
— pl; | pas 1º + sly — gl; 

1,0 0,0208 0,1146 0,0208 0,0296 0,0573 0,0088 0,0088 0,0573 0,0296 
1,1 0,0189 0,1272 0,0270 0,0301 0,0546 0,0084 0,0112 0,0726 0,0354 
Tod 0,0166 0,1421 0,0335 0,0304 0,0521 0,0080 0,0138 0,0900 0,0415 
1.3 0,0140 0,1593 0,0405 0,0308 0,0498 0,0077 0,0168 0,1095 0,0482 
1,4 | 0,0110 0,1787 0,0479 0,0311 0,0477 0,0073 0,0202 0,1310 0,0552 
1,5 0,0076 0,2005 0,0557 0,0314 0,0458 0,0071 0,0238 0,1547 0,0627 
1,6 | 0,0039 0,2246 0,0639 0,0317 0,0441 0,0068 0,0278 0,1805 | 0,0707 
1,7 — 0,0001 | 0,2509 0,0725 0,0319 0,0424 0,0065 0,0321 0,2085 0,0791 
1,8 — 0,0046 | 0,2796 0,0816 | 0,0322 0,0409 0,0063 0,0367 0,2387 0,0879 
1,9 — 0,0093 | 0,3105 0,0911 0,0324 0,0395 0,0061 0,0417 0,2710 0,0972 
2,0 — 0,0144 | 0,3437 0,1010 0,0326 0,0382 0,0059 0,0470 0,3056 0,1068 
2,25 — 0,0287 | 0,4368 0,1275 0,0330 0,0353 0,0054. 0,0618 0,4016 0,1329 
2,5 — 0,0465 | 0,5443 0,1567 0,0334 0,0327 0,0050 0,5115 0,1617 


la/ 1, 

1,0 0,0417 
11 0,0384 
Ta 0,0345 
1,3 | 0,0301 
1,4 | 0,0250 
1,5 0,0193 
1,6 0,0131 
1,7 0,0063 
1,8 — 0,0012 
1,9 — 0,0092 
2,0 — 0,0179 
EE A — 0,0420 


2,5 | — 0,0699 
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0,1042 
0,1156 
0,1292 
0,1448 
0,1625 
0,1823 
0,2042 
0,2281 

0,2542 
0,2823 
0,9125 
0,3971 


0,4948 


0,0417 
0,0534 
0,0660 
0,0793 
0,0936 
0,1086 
0,1245 
0,1412 
0,1588 
0,1772 
0,1964 
0,2481 
0,3051 


0,0565 
0,0573 
0,0579 
0,0585 
0,0590 
0,0595 
0,0600 
0,0604 
0,0608 


0,0612 


0,0615 
0,0623 


0,3136 


0,0521 
0,0496 
0,0473 
0,0453 
0,0434 
0,0417 
0,0401 
0,0386 
0,0372 
0,0359 
0,0347 
0,0321 
0,0298 


0,0149 
0,0142 
0,0135 
0,0129 
0,0124 
0,0119 
0,0114 
0,0110 
0,0106 
0,0103 
0,0099 
0,0092 
0,0085 


| 0,0787 


0,0149 
0,0189 
0,0234 
0,0284 
0,0340 
0,0402 
0,0469 
0,0542 
0,0620 
0,0704 
0,0794 
0,1043 


0,1329 


0,0521 
0,0660 
0,0818 
0,0995 
0,1191 
0,1406 
0,1641 
0,1895 
0,2170 
0,2464 
0,2778 
0,3651 


0,4650 


0,0565 
0,0676 
0,0795 
0,0923 
0,1060 
0,1205 
0,1360 
0,1523 
0,1694 
0,1875 
0,2063 
0,2573 
0,0629 
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TABELA 6 — Viga continua de dois tramos (cont.) 


n=3/4 


bh EE ella ig! 
= lo + sly = "Bla: 
O 0,0625 0,0938 0,0625 0,0813 0,0469 0,0188 0.0188 0,0469 0,0813 
1,1 0,0584 0,1041 0,0794 0,0821 0,0446 0,0179 0,0238 0,0594 0,0972 
1,2 0,0535 0,1163 0,0975 0,0830 0,0426 0,0170 0,0295 0,0736 0,1145 
Td 0,0479 0,1303 0,1169 | 0,0837 0,0408 0,0163 0,0358 0,0896 0,1332 
1,4 0,0415 0,1463 0,1375 | 0,0844 0,0391 0,0156 0,0429 0,1072 0,1531 
1,8 0,0344 0,1641 0,1594 0,0850 0,0375 0,0150 0,0506 0,1266 0,1744 
1,6 0,0265 0,1838 0,1825 0,0856 0,0361 0,0144 0,0591 0,1477 0,1969 
1! 0,0179 0,2053 0,2069 0,0861 0,0347 0,0139 0,0682 0,1706 0,2208 
1,8 0,0085 0,2288 0,2325 0,0866 0,0335 0,0134 0,0781 0,1953 0,2459 
1,9 0,0016 | 0,2541 0,2594. 0,0871 0,0323 0,0129 0,0887 0,2217 0,2723 
2,9 — 0,0125 | 0,2813 0,2875 0,0875 0,0312 0,0125 0,1000 0,2500 0,3000 
2,25 — 0,0430 | 0,3574 0,3633 0,0885 0,0288 0,0115 0,1314 0,3286 0,3748 
2,5 — 0,0781| 0,4453 0,4409 0,0893 0,0268 0,0107 0,1674 0,4185 0,4576 
n=1 
|:/1, kt | k k! k k k' ko k 
— plé = 8h, + sly = 
LO 0,0833 0,0833 0,0833 0,1042 0,0417 0,0208 0,0208 0,0417 0,1042 
1,1 0,0787 0,0925 0,1050 0,1052 0,0397 0,0198 0,0264 0,0528 0,1248 
Tiniê 0,0733 0,1033 0,1283 0,1061 0,0379 0,0189 0,0327 0,0655 0,1473 
1,3 0,0671 0,1158 0,1533 0,1069 0,0362 0,0181 0,0398 0,0796 0,1714 
1,4 0,0600 0,1300 0,1800 0,1076 0.0347 0,0174 0,0476 0,0953 0,1974 
1,5 0,0521 0,1458 0,2083 0,1083 0,0333 0,0167 0,0562 0,1125 0,2250 
1,6 0,0433 0,1633 0,2383 0,1090 0,0321 0,0160 0,0656 0,1313 0,2544 
1,7 0,0337 0,1825 0,2700 0,1096 0,0309 0,0154 0,0758 0,1516 0,2854 
1,8 0,0233 0,2033 0,3033 0,1101 0,0298 0,0149 0,0868 0,1736 0,3182 
1,8 0,0121 0,2258 0,3383 0,1106 0,0287 0,0144 0,0985 0,1971 0,3527 
2,0 0,0000 0,02500 0,3750 0,1111 0,0278 0,0139 0,1111 0,2222 0,3889 
2,25 -— 0,09339 | 0,3177 0,4740 0,1122 0,0256 0,0128 0,1460 0,2921 0,4868 
0,1131 0,0238 0,0119 0,1860 0,3721 
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2,5 


|— 0,0729 | 0,3958 


0,5833 


TECNICA 424 


0,5952 


TABELA 7 — Viga continua de três tramos parcialmente encastrada nos apoios extremos, com cargas uniforme- 
mente distribuidas 


(n) p+s (n) 
7 VOO RNOAADARAOSAROANAGDANAAAAAGHAGARACEARARAARAARANAA 4 
A AB Ye: SD 


pmiguece lo e 


lo/1, 


— pl? 
0,4 0,0000 0,0831 0,0000 0,0911 | 0,0130 0,0000 0,0000 0,0050 
0,5 0,0000 0,0804 0,0000 0,0857 0,0143 0,0000 0,0000 0,0089 
0,6 0,0000 0,0800 0,0009 0,0810 0,0152 0,0000 0,0000 0,0142 
0,7 0,0000 0,0819 0,0000 0,0768 0,0158 0,0009 0,0000 0,0209 
0,8 0,0000 0,0859 0,0000 0,0731 0,0162 0,0000 0,0000 0,0291 
0,9 0,0000 0,0929 0,0000 0,0697 0,0165 0,0009 0,0000 0,0388 

1,0 0,0000 0,1000 0,0009 0,0667 0,0167 0,0000 0,0000 0,0500 
1,1 0,0000 0,1100 0,0000 0,0639 0,0167 0,0000 0,0000 0,0628 
1,2 0,0000 0,1218 0,0000 0,0614 0,0167 0,0000 0,0000 0,0771 
1,25 0,0000 0,1284 0,0000 0,0602 0,0167 0,0000 0,0000 0,0849 
1,3 0,0000 0,1355 0,0000 0,0591 0,0167 0,0000 0,0000 0,0931 
1,4 0,0000 0,1510 0,0000 0,0569 0,0166 0,0000 0,0000 0,1106 
1,5 0,0000 0,1683 J,0000 0,0549 0,0165 0,0000 0,0000 0,1298 
1,6 0,0000 0,1874 0,0000 0,0531 0,0163 0,0000 0,0000 0,1506 
17 0,0000 0,2082 0,0900 0,0514 0,0162 0,0000 0,0000 0,1730 
1,8 0,0000 0,2308 0,0000 0,0498 0,0160 0,00900 0,0000 0,1970 
1,9 0,0000 0,2552 0,0000 0,0483 0,0158 0,0000 0,0000 0,2227 
2,0 0,0000 0,2812 0,0000 0,0469 0,0156 0,0000 0,0000 0,2500 
2,25 0,0000 0,3540 0,0000 0,0437 0,0151 0,0000 0,0000 0,3254 
2,5 0,0000 0,4375 0,0000 0,0409 0,0146 0,0000 0,0000 0,4112 
3,0 0,0000 0,6364 0,0000 0,0364 0,0136 0,0000 0,0000 0,6136 
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TE sl,” 


sl, 


1 
= s| * 
| 


TABELA 7 — Viga continua de três tramos (cont. 


p+*s 


204 


1/1, p= E. | ka=k, | ki=ky | li=k =k" | ko=hk, lee | et = 
= pl — sl, sl; — sly sh — ly 
0,4 0,0270 0,0747 0,0259 0,0818 | 0,0124 0,0019 0,0008 0,0053 
0,5 0,0273 0,0726 0,0267 0,0767 0,0135 0,0021 0,0014 0,0093 
0,6 0,0273 0,0728 0,0273 0,0722 0,0142 0,0022 0,0023 0,0148 
0,7 0,0269 0,0753 0,0279 0,0683 0,0147 0,0023 0,0033 0,0217 
0,8 0,0262 0,0800 0,0285 0,0649 0,0151 0,0023 0,0046 0,0301 
0,9 0,0251. 0,0866 0,0290 0,0618 0,0152 0,0023 0,0062 0,0401 
1,0 0,0238 | 0,0952 0,0294 0,0590 0,0153 0,0024 0,0079 0,0516 
1,1 0,0222 0,1058 0,0298 0,0565 0,0153 0,0024 0,0099 0,0647 
1a 0,0203 0,1182 0,0301 0,0541 0,0153 0,0024 0,0122 0,0793 
1,25 0,0192 0,1251 0,0303 0,0531 0,0153 0,0023 0,0134 0,0873 
1,3 0,0181 0,1324 0,0305 0,0520 0,0152 0,0023 0,0147 0,0956 
1,4 0,0156 | 0,1484 0,0308 0,0501 0,0151 0,0023 0,0175 0,1135 
1,5 0,0129 0,1663 0,0310 0,0483 0,0149 0,0023 0,0205 0,1330 
1,6 0,0099 0,1859 0,0313 0,0466 0,0148 0,0023 0,0237 0,1541 
Tok 0,0066 0,2073 0,0315 0,0451 0,0146 0,0022 0,0272 0,1768 
1,8 0,0030 0,2304 0,0318 0,0436 0,0144 0,0022 0,0310 0,2012 
13 — 0,0008 0,2553 0,0320 0,0423 0,0142 0,0022 0,0350 0,2272 
2,0 — 0,0049 0,2819 0,0322 0,0410 0,0149 0,0022 0,0392 0,2549 
2,25 — 0,0163 0,3559 0,0326 0,0381 0,0135 09,0021 0,0510 0,3313 
Lo — Q,UZIS 0,4406 0,0339 0,0356 0,0130 0,0020 0,0643 0,4180 
3,0 — 0,0603 0,6418 0,0336 0,0316 0,0121 0,0019 0,0957 0,6223 


A puto e | 
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Á 


TABELA 7 


a/14 a=k, | k,=k 
— pl? 
0,4 0,0524 | 0,0668 
0,5 0,0528 | 0,0653 
0,6 0,0525 | 0,0662 
0,7 0,0516 | 0,0693 
0,8 0,0501 | 0,0745 
0,9 0,0481 | 0,0817 
1,0 0,0455 | 0,0909 
Li 0,0423 | 0,1020 
1,2 0,0386 | 0,1149 
1,25 0,0366 | 0,1220 
1,3 0,0344 | 0,1296 
1,4 0,0297 | 0,1462 
1,5 0,0244 | 0,1645 
1,6 0,0187 | 0,1846 
1,7 0,0124 | 0,2065 
1,8 0,0057 | 0,2300 
1,9 — 0,0015 | 0,2554 
2,0 — 0,0093 | 0,2824 
2,25 — 0,0307 | 0,3575 
2,5 pe 0,0552 | 0,4433 
3,0 — 0,1133 | 0,6467 
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ki=k! 


0,0506 
0,0520 
0,05391 
0,0542 
0,0551 
0,0559 
0,0566 
0,0573 
0,0579 
0,0582 
0,0584 
0,0589 
0,0594 
0,0598 
0,0602 
0,0606 
0,0609 
0,0613 
0,0620 
0,0626 


0,0636 


= s);º 


ot = ke 


[E E 


0,0729 


0,0681 


0,0640 


0,0604 


0,0572 


0,0544 


0,0518 


0,0495 


0,0474 


0,0465 


0,0455 


0,0438 


0,0421 


0,0406 


0,0393 


0,0380 


0,0367 


0,0356 


0,0331 


0,0309 


0,0273 


0,0116 
0,0125 
0,0132 
0,0136 
0,0138 
0,0139 
0,0140 
0,0139 
0,0138 
0,0138 
0,0137 
0,0136 
0,0134 
0,0132 
0,0130 
0,0129 
0,0127 
0,0125 
0,0120 
0,0115 


0,0106 


— Viga continua de três tramos (cont.) 


(n) 


MOQURAGERANAGAO OO OENGRAAANGNAGANARANANAO, c 


— sly 
0,0033 
0,003 
0,0038 
0,0039 
0,0039 
0,0040 
0,0040 
0,0040 
0,0040 
0,0039 
0,0039 
0,0039 
0,0038 
0,0038 
0,0037 
0,0037 
0,0036 
0,0036 
0,0034 
0,0033 


0,0030 


ké =! 
sly 
0,0016 
0,0028 
0,0044 
0,0064 
0,0089 
0,0118 
0,0152 
0,0190 
0,0232 
0,0255 
0,0280 
0,0331 
0,0388 
0,0449 
0,0515 
0,0586 
0,0661 
0,0741 
0,0961 
0,1211 


0,1800 


kg = ka 


— sly 
0,0055 
0,0097 
0,0154 
0,0225 
0,0311 
0,0413 
0,0530 
0,0664 
0,0813 
0,0894 
0,0978 
0,1160 
0,1358 
0,1572 
da BO2 
0,2049 
0,2313 
0,2593 
03364 
0,4239 


0,6300 
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to /1, 


2,25 


TABELA 7 — Viga continua de três tramos (cont.) 


(n) 


Á 
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n=3/4 
=k, [kt | M=6p [RSA g | R=kpP | SD | ki=ho | kiko 
— pl — 81º sl, E — sl; sl; | — gl; 

0,0763 0,0593 0,0743 0,06 0,0107 0,0043 0,0023 0,0057 
0,0766 | 0,0585 0,0760 0,0599 0,0115 0,0046 0,0040 0,0101 
0,0760 | 0,0600 0,0775 0,0561 0,0120 0,0048 0,0064. 0,0159 
0,0745 0,0637 0,0788 0,0529 0,0123 0,0049 0,0093 0,0232 
0,0722 0,0695 0,0800 0,0500 0,0125 0,0050 0,0128 0,0320 
0,0691 0,0773 0,0810 0,0474 0,0126 0,0050 0,0170 0,0424 
0,0652 | 0,0870 0,0819 0,0452 0,0125 0,0050 0,0217 0,0543 
0,0606 0,0985 0,0828 0,0431 0,0125 0,0050 0,0272 0,0679 
0,0552 0,1119 0,0835 0,0412 0,0124 0,0049 0,0332 0,0831 
0,0523 0,1193 0,0839 0,0403 0,0123 0,0049 0,0365 0,0913 
0,0491 0,1271 0,0842 0,0395 0,0122 0,0049 0,0399 0,0999 
0,0423 0,1441 0,0848 0,0379 0,0121 0,0048 0,0473 0,1183 
0,0348 0,1629 0,0854 0,0365 0,0119 0,0048 0,0553 0,1383 
0,0266 | 0,01834 0,0859 0,0352 0,0117 0,0047 0,0640 0,1600 
0,0177 0,2057 0,0864 0,0339 0,0115 0,0046 0,0733 0,1833 
0,0081 | 0,2297 0,0869 0,0328 0,0113 0,0045 0,0833 0,2083 
— 0,0022 | 0,2554 00873 0,0317 0,0112 00045 0,0940 0,2349 
— 0,0132 | 0,2829 0,0877 0,0307 0,0110 0,0044 0,1053 0,2632 
— 0,0436 0,3590 0,0886 0,0285 0,0105 0,0042 0,1364 0,3410 
| — 0,0783 0,4457 0,0894 0,0265 0,01019 0,0040 0,1717 0,4293 
— 0,1604 0,6509 0,0906 0,0234 0,0092 0,0037 0,2547 0,6363 


TABELA 7 — Viga contínua de três tramos (cont.) 


(n) p+s (n) 


A 
A! A B AG AS 


q — 


| 


A 


n= 1 

71, | to | =% | ki=kr | ki=ky | k=ky ki =k! 

— ply iii sl, -— dy sl; -— al; 
0,4 0,0989 0,0552 0,0970 | 0,0560 0,0097 0,0049 0,0030 0,0059 
0,5 0,0990 | 0,0521 0,0990 0,0521 0,0104 0,0052 0,0052 0,0104 
0,6 0,0979 | 0,0542 0,1006 0,0487 0,0108 0,0054 0,0082 0,0164 
0,7 0,0957 | 0,0585 01021 0,0458 0,0110 0,0055 0,0119 0,0238 
0,8 0,0926 | 0,0649 0,1034 0,0432 0,0111 0,0056 0,0164 0,0328 
0,9 0,0884 | 0,0732 0,1045 0,0409 0,0112 0,0056 0,0217 0,0434 
1,0 0,0833 | 0,0833 0,1056 0,0389 0,0111 0,0056 0,0278 0,0556 
1,1 0,0773 | 0,0954 0,1065 0,0371 0,0110 0,0055 0,0347 0,0693 
1,2 00704 | 0,1092 0,1073 0,0354 0,0109 0,0054 0,0424 0,0847 
1,25 0,0666 | 0,1168 0,1077 0,0346 0,0108 0,0054 0,0465 0,0930 
1,3 0,0626 | 0,1249 0,1081 0,0339 0,0107 0,0054 0,0509 0,1017 
1,4 0,0539 | 0,1423 0,1087 0,0325 0,0106 0,0053 0,0602 0,1204 
1,5 0,0443 | 0,1615 0,1094 0,0312 0,0104 0,0052 0,0703 0,1406 
1,6 0,0338 | 0,1824 0,1100 0,0301 0,0102 0,0051 0,0813 0,1625 
13 0,0225 | 0,2050 0,1105 0,0290 0,0101 0,0050 0,0930 0,1861 
1,8 0,0103 | 0,2294 0,1110 0,0280 0,0099 0,0049 0,1057 02113 
1,9 | — 0,0028 | 0,2555 0,1115 0,0271 0,0097 0,0049 0,1191 0,2382 
2,0 — 0,0167 | 0,2833 0,1119 0,0262 0,0095 0,0048 0,1333 0,2667 
2,25 — 0,0552 | 0,3603 0,1129 0,0242 0,0091 0,0045 0,1726 0,3452 
2,5 — 0,0990 | 0,4479 0,1137 0,0226 0,0087 0,0043 0,2170 0,4340 
3,0 — 0,2024 | 0,6548 | 0,1151 0,0198 0,0079 0,0040 0,3214 0,6429 
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LUOS 


0000'0 
0000'0 
0000'0 
0000 0 
0000'0 
0000'0 
0000'0 
0000'0 
0000'0 
0000'0 
0000'0 
00000 
00000 
0000'0 
00000 
0000 0 
0000'0 
0000'0 
0000'0 
0000'0 


0000'0 


p= 


a 


“SOLUBI]XO 


EZ9LO 
Z60L O 
99800 
L990'0 
v6S0 0 
9250 0 
c9t0 O 
EOrO O 
LvECGO 
9620'0 
6rZo O 
BLTO O 
LOZO'O 


89L0'0 


vELOO 
POLC'O 
8.000 
9500 0 
8E00 0 
FZOo'o 


ELCOO 


: 18 + 


solode sou 


LTEP O 
LSCE O 
ECST o 
0002'0 
PLBLO 
LE9LC 
B94L'O 
BOELO 
BSLLO'O 
9LOLO 
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Espectros de massa e energia na formação de 


películas finas usando o processo de 


"evaporação por impacto” 


RESUMO 


O conhecimento dos espectros de massa e de 
energia das partículas incidentes no substrato, quando 
da formação de Películas Finais é de extrema impor- 
tância para a compreensão dos seus processos de for- 
mação e crescimento, intimamente ligados às proprie- 
dades macroscópias que apresentam. Para a determi- 
nação destes espectros, construiu-se um analisador de 
massa e de energia cuja montagem se adaptou ao estudo 
do processo de «Evaporação por Impacton. 


1 — INTRODUÇÃO 


«Pliculas finas» são camadas de material deposi- 
tado sobre uma base adequada, ou substrato, com espes- 
suras variando de poucas centenas de À” a alguns mi- 
crons, o que corresponde a dimensões da ordem de 
algumas camadas atómicas. 


O uso de «Películas Finas» tem sido largamente 
conhecido desde há centenas de anos, mas só muito 
recentemente se reconheceu o seu interesse em domi- 
nios que englobam quer a obtenção de sistemas ópticos 
e magnéticos, quer a electrónica duma maneira geral 
e a microslectrónica duma forma especial, quer ainda 
o estudo das estruturas de sólidos e suas propriedades 
fisicas, naturalmente acompanhando o progresso e exi- 
gências inerentes âqueles domínios. 


Estas películas, são produzidas depositando mate- 
rial sobre um suporte (base ou subtrato) utilizando-se 
para isso diferentes técnicas, algumas das quais actual- 
mente em facto ds pleno desenvolvimento. O crescimento 
da película durante a deposição envolve usualmente um 
estado de nucleação que caracteriza a formação de grupos 
estáveis de átomos e a subsequente coalescência ou 
união desses grupos. Às características finais da peli- 
cula, serão função de vários parâmetros, entre os quais 
se destacam a temperatura do substrato, a velocidade 
de deposição e determinados parâmetros energéticos. 


Existem várias razões para a constante expansão 
do uso de «Películas Finas» na Microelectrónica, das 
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ABSTRACT 


The nature of the particles incident on a substracte 
is od basic interest for on understanding of the nuclea- 
tion process of thin films obtained by deposition. 

A sintable device has been constructed to study 
the nature and energy of those particles, using a «catho- 
dic Sputtering» process. 


quais se evidência a possibilidade que estas películas 
oferecem em se conseguir circuitos de pequenas dimen- 
sões, a facilidade de produção de substâncias de elevada 
pureza ou de materiais de composição rigorosamente 
controlada e ainda a possibilidade de obtenção de 
circuitos com determinadas características de funda- 
mental importância. No que respeita a tentativas já 
efectuadas para a utilização de Películas Finas como 
únicos elementos formadores de circuitos, tendo em 
vista as vantagens daí resultantes ao baixo custo de 
fabricação e facilidade de materialização de esquemas 
dificeis, não tem sido possível até hoje usar desse 
benefício, pois não se conseguiu ainda suplantar a 
actual tecnologia de circuitos monolíticos. No entanto, 
interessantes desenvolvimentos parecem possíveis sendo 
licito esperar rápidos progressos, pois a citada inibição 
resulta principalmente do atraso relativo em que a 
técnica de crescimento de Películas Finas se encontra 
actualmente e não possibilidades que à primeira vista 
O processo em si oferece, 


2 — À «EVAPORAÇÃO POR IMPACTO» NA FORMAÇÃO 
DE PELÍCULAS FINAS 


Este método, traduz o fenómeno da ejecção de par- 
ticulas de energia suficiente, O material libertado, na 
forma atómica ou molecular, pode condensar-se sobre 
uma superfície, formando depósitos com as propriedades 
características das Películas Finas obtidas por outros 
meios, 
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A realização laboratorial do processo, requer O 
bombardeamento dum alvo por partículas com suficiente 
energia, as quais podem ser particulas neutras ou iões. 
A simplicidade em acelerar partículas ionizadas e a 
dificuldade em obter particulas neutras com energia 
cinética suficiente, levaram ao uso sistemático de iões 
positivos, os quais podem ser gerados por dois processos 
fundamentais. 


(1) Através duma descarga luminosa em meio 
gasoso ionizado, como o alvo servindo simul- 
tâneamente de eléctrodo negativo, podendo 
o eléctrodo colector coincidir ou não com 
o ânodo. Neste caso, é aplicada uma tensão 
continua de valor suficiente para desenca- 
dear e manter o fenómeno; 


(2) — Usando um sistema que produza um feixe 
de iões que bombardeiam directamente o 
alvo. 


Poce utilizar-se por exemplo uma cor- 
rente de electrões provenientes dum cátodo 
quente, a qual passando por um meio 
gasoso, origina a formação dum plasma. 


Os iões são extraidos e acelerados na di- 
recção do alvo, por acção dum campo eléc- 
trico independente. 


Extractor 
de 
particulas 


Cámara de 


Anoalizador 
de 
energia 


| evaporação 


por impacto 


l 
I 
I 
I 
| 
I 
I 
L. 


Fig. 1 — Diagrama de blocos 


No que respeito propriamente à deposição, ou seja 
a nucleação e crescimento da película, foi constatado 
que ocorre diferentemente de acordo com o processo 
utilizado, De facto, tomando os dados publicados recen- 
temente por vários autores , o uso da evaporação por 
impacto através duma descarga luminosa em meio 
gasoso ionizado, implica na primeira fase de cresci- 
mento das películas uma densidade de ilhas superior 
a encontrada para o processo de evaporação térmica 
no vácuo e uma maior facilidade de crescimento epi- 
taxial, não tendo sido no entanto dada explicação 
plausível para esta ocorrência, Entre os factores que 
poderão contribuir para esta diferença de comporta- 
mento, tomam vulto a natureza e a energia que caracte- 
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rizam as peliculas incidentes numa pelicula na fase de 
formação, À comprovar esse facto, pode indicar-se os 
trabalhos de Layton e Campbel(') donde se conclui 
que o crescimento epitaxial das películas depositadas 
está associado com a energia das partículas incidentes 
e o trabalho de Stirland(') que mostrou que o bombar- 
deamento electrónico do substracto pode desempenhar 
um papel importante no crescimento da película fina. 
Tem também sido sugerido que átomos ou electrões 
incidentes, podem produzir defeitos pontuais do subs- 
tracto com resultado da sua penetração, podendo esses 
de'eitos actuar como pontos de partida para a formação 
de núcleos estáveis, 


3 — DESCRIÇÃO DO TRABALHO A REALIZAR 


Tomando como que uma informação detalhada 
acerca da massa e energia das películas incidentes na 
superficia onde a pelicula cresce, poderá explicar as 
razões que determinam as características particulares 
de cada particula depositada, projectâmos e construímos 
um sistema de modo a permitir essa observação('). O 
sistema apresenta-se em diagrama na figura 1, seguindo 
em certa medida o processo usado por J. Coburnl') que 
apresentou pela primeira vez um trabalho especifica- 
mente devotado à determinação do espectro de massa 
relativa ao fenómeno da «Evaporação por Impacto» 


Acelzrador Analizador | Collector de | 
iões e 
aparelhagem 


de medida 


de | de 
particulas 


do analizador de massa e energia 


usando uma descarga em meio gasoso ionizado, com 
tensão continua aplicada aos electrodos. 


A primeira parte do trabalho incidiu apenas sobre 
a determinação do espectro de energia das partículas 
incidentes e dentro destas somente se detectaram 
as particulas carregadas, eventualmente iões e electrões, 
pois a observação de partículas neutras exigira a intro- 
dução duma câmara de ionização por impacto, a qual 
se planeia instalar numa fase mais adiantada do projecto. 
O processo de deposição usado foi o de Evaporação 
por Impacto, utilizando-se o ÁArgon como gás ambiente, 
com tensão continua aplicada aos electrodos. 

A figura 2 representa esquematicamente a dispo- 
sição seguida. 
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entrada de gás 


&nalisador 
de energia 


substrato & 
lextractor de 
partículas 


câmara de pulverização catódica 


| câmara de vácuo 


bomba de difusão de óleo 


Fig 2 — Representação esquemática do arranjo para a deter- 
minação do espectro de energia das partículas incidentes no 
substrato durante a formação duma pelicula fina 


Seguindo a análise de Purcell(*), obtém-se um 
efeito de focagem das orbitas da partículas que saindo 
dum dado ponto (fonte) e entrada no condensador 
segundo um feixe divergente, são agrupadas após uma 
rotação de 180º de acordo com a energia que possuem 
e para para dada diferença de potencial V, aplicada às 
armaduras do condensador. De acordo com a figura 3, 
tem-se: 


onde E representa a energia da particula (em volts 
equivalentes, qE — electrão volts) que segue a direcção 
normal na região (7) 


O analisador de energia utilizado, é do tipo elec- 
trostático, de deflexão, consistindo fundamentalmente 
num sector de 90º formado por duas superfícies esfé- 
ricas concêntricas. 


Segue-se portanto que todas as partículas duma 
dada energia E deixando A segundo um feixe divergente 
de pequeno a, serão aproximadamente focadas em B, 
desde que se aplique o valor correcto da diferença de 
potencial V, entre as duas superfícies esféricas. 
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Fio. 3 —- Representação bidimensional do analisador de energia 


Construiu-se o sistema, usando dois sectores esfé- 
ricos em Alumínio com & =— 90º e de diâmetros res- 
pectivamente 10 cm e 11 cm, donde resulta a = b = ,, 2 
Fr, sendo r; = 5,25 cm, Estas dimensões conduzem à 
relação teórica entre a tensão de focagem e a energia 
das partículas: 


Ve 
o = 104 
E 


A figura 4 representa por sua vez um pormenor do 
extractor de partículas que se incluiu na parte inferior 
da câmara de deposição, entre esta e o analisador de 
energia. 


5 


1- Orificio de extracção 
2- Basc ou substrato 


3- Tampa de topo da câmara de pulverização cotódica 
4- Analizador de energia 
5- Anel isolante 


Fig. 4 — Extractor de partículas 


O elemento que contém o orifício de extracção, 
está eléctricamente isolado, quer do analisador de 
energia, quer do substracto, de modo a permitir a obser- 
vação do efeito provocado pela variação da sua polariza- 
ção relativamente ao resto do substrato, o qual funciona 
ao mesmo potencial do ânodo. De facto, sendo aquele 
elemento concêntrico com o substrato e ocupando ape- 
nas uma pequena parte da área deste, funcionará ele 
próprio como substrato, sempre que o seu potencial é 
diferente do potencial do ânodo. 


O colector de particulas consiste simplesmente 
num alvo de aluminio de forma cilíndrica, com uma 
pequena abertura para entrada das partículas e directa- 
mente ligado ao circuito dum electróômetro. À fim de 
evitar os efeitos da emissão secundária, incluiu-se tam- 
bém um supressor ligado a um potencial negativo de 


215 


valor pequeno, de modo a não provocar deflexão dos 
iões primários, 

Duma maneira resumida, o funcionamento do 
sistema, processa-se da seguinte forma: um feixe de 
partículas provenientes da descarga, é extraído através 
do orifício aberto no substrato, entrando no analisador 
de energia instalado numa câmara onde a pressão se faz 
suficientemente baixa, As partículas separadas no ana- 
lisador de energia de acordo com o valor energético 
pretendido, embatem directamente no alvo do colector 
de particulas, originando uma corrente, a qual é lida 
através dum electrómetro, 


4 -— RESULTADOS EXPERIMENTAIS 


Durante o primeiro grupo de experiência, usou-se 
um orifício de extracção de partículas, de 1,8 mm de 
diâmetro o que permitiu um ângulo de aceitação de 6,5.º 
Argon com gás ambiente na câmara de Evaporação por 
Impacto, variando-se a sua pressão de 404 a 9Du 
(1u = 10" mmHg), o que levou uma pressão na câmara 
do analisador de energia de aproximadamente 10 “STorr, 
valor considerado suficiente para as experiências levadas 
a cabo. À tensão entre eléctrodos variou-se de 600 
volts a 1.500 volts. 

O material escolhido para o alvo foi o cobre. À gama 
de utilização da tensão de Evaporação por Impacto 
exige que a tensão de focagem a aplicar ao analisador 
de energia, se situe dentro dos limites: 

Vp=o volts, V,5= 0,191 X 1500 volts 
de modo a que se possa varrer todos os valores possi- 
veis da energia das particulas que chegam ao substrato. 


A tensão a aplicar aos sectores esféricos interior e exte- 
rior será: 


3,191 
Ve Volts 


respectivamente ce valor positivo e negativo no caso 
de se pretender sintonizar partículas negativas e de 
polaridade contrária se se pretender partículas positivas, 
isto para se conservar sempre a parte central do anali- 
sador ao potencial zero, 


As curvas da figura 5 mostram a distribuição ener- 
gética das particulas incidentes no substrato, para uma 
primeira séria de experiências, cujos parâmetros são 
os seguintes: 


- Analisador de energia sintonizado sempre para 
partículas negativas, 


(a) Angulo de aceitação y = 6,5 
Tensão aplicada aos eléctrodos V — 1500 


volts 
Corrente total de descarga | = 13 mA 
Pressão do gás ambiente p= 53u 
Tensão de polarização do substrato Va — O 
volts 
(bh) 
y = 6,5" 
V = 1.090 volts 
|, = 6 mÃ 
p — 53u 
VB = q volts 
(c) 
Y = 6,5º 
V = 1.000 volts 
ly = 11 mA 
p = 70u 


Wa Õ volts 


E sui | 7 
Eta zona de baixos energios Ts zona de elevadas energias 
nerg T T | q Ea =1205V] |V=1445V 
S p=45u p= 45 
28 p= 45 
26 
24 
ess Ze 
20 20 
18 om 18 
E 3 16 
4 y 14 
14454 
12 dp R p= op 12 
[1 
to 1 ã 10 
] V=1205V ú 
OB dm p=45u A + 
ls 6! 
beba da | 
CF 1 DE a | 
04 cd d 
j 
02 i 2 
E-=Ve/O1M eV 
—A+— ++ dy —+— |-—— 4 = 
O 2 4 6 B 10 20 30 50 Cc 120 200 300 Ve lvolts) 
Fig. 5 — Curvas de distribuição de energia 
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Para se facilitar a representação gráfica, desta- 
caram-se duas regiões correspondentes respectivamente 
à zona das baixas energias e à zona das elevadas ener- 
gias, que na realidade correspondem exactamente às 
duas contribuições importantes dentro de toda a variação, 
conforme experimentalmente se observou, Em ordenadas 
fez-se a representação da corrente normalizada, isto é 
valores da corrente medida através do electrómetro 
divididos pela corrente total da descarga, usando-se duas 
escalas diferentes conforme se pode constatar na figura. 


5 — INTERPRETAÇÃO DOS RESULTADOS 


Da observação das curvas obtidas, pode concluir-se 
em princípio a existência das seguintes tendências: 


a) A distribuição em energia apresenta essencial- 
mente duas contribuições distintas: a relativa à região 
de baixas energias, com máximos na zona dos 8a 12 eV 
estendendo-se numa gama de alguns eV e a relativa à 
região de elevadas energias, com valores da corrente 
crescendo até ao máximo valor da tensão aplicada (em 
equivalentes eV) onde se observa práticamente um 
corte, qualquer que seja essa tensão aplicada, sendo o 
corte mais pronunciado à medida que a tensão aplicada 
aumenta de valor. Tudo indica que nesta região existe 
um limiar de energia, que determina a existência de 
particulas chegando ao substracto. 


bj) À posição dos máximos na zona de baixa energia 
não varia pronunciadamente quer com a tensão apli- 
cada, quer com a pressão do gás ambiente. 


c) 4 amplitude dos máximos, por sua vez, depende 
Gecididamente da pressão reinante, mais pronuncia- 
damente que da tensão aplicada e a extensão em ener- 
gia coberta pelo ramo descendente da curva, é maior 
para pressões mais elevadas. 


d) O corte observado na região de elevadas ener- 
gias, depende unicamente da tensão aplicada, no que 
respeita à sua posição relativamente à energia, mas a 
amplitude das respectivas curvas, depende efectivamente 
da pressão do gás ambiente, 


e) Quanto mais se aumenta a tensão aplicada, mais 
os valores achados se concentram junto dos respectivos 
máximos das curvas, tornando-se o máximo cada vez 
mais pronunciado pelo lado ascendente das curvas na 
direcção das energias crescentes. Tudo indica, que nes- 
tas condições, o feixe de particulas tende a tornar-se 
monocenergético, de valor energético igual ao corres- 
pondente à tensão aplicada, em equivalentes ev, 


Tomando em consideração quer a teoria geral das 
descargas luminosas em meios gasosos,quer os resul- 
tados de determinados autores, conseguidos através 
de experiências elaboradas em diferentes condições de 
trabalho, quer ainda a nossa própria experiência, criâmos 
o seguinte modelo para interpretação dos resultados 
obtidos: 
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a) O máximo que se situa na parte extrema da re- 
gião de elevadas energias, representa os electrões que 
formando um feixe quase monoenergético, atravessam 
toda a distância que separa o cátodo do ânodo sem 
sofrerem perdas apreciáveis devidas a colisões. Estes 
electrões são formados imediatamente junto do cátodo. 


b) O espectro de energia intermédio, representa 
os electrões que tendo sido criados através da descarga, 
sofrerem bastantes colisões não elásticas, 


c) O máximo situado na região de pequenas ener- 
gias, representa os electrões que emergindo da coluna 
positiva animados duma pequena energia cinética, são 
atraídos para o ânodo através do espaço escuro anódico. 


Este máximo poderá também ser devido a iões 
negativos do gás ambiente, ou a iões negativos do metal 
do alvo (cátodo) ejectados deste como partículas neu- 
tras, mas tendo sofrido uma posterior ionização ao atra- 
vessarem o meio gasoso ambiente, (ionização por 
impacto). 


A primeira destas hipóteses referentes a iões, é 
no entanto bastante improvável, porque o gás utilizado 
(Argon) tem uma área de captura para ionização muito 
pequena (a área de captura para ionização é definida 
como o inverso do percurso médio livre). 

A segunda das hipóteses poderá todavia ser consi- 
derada como válida 


A fim de justificar o modelo apresentado, torne- 
mos presente o que acontece no espaço escuro catódico 
numa descarga luminosa «Self-mantidas. Os electrões 
que atravessam aquela região, são constituídos por dois 
grupos distintos, ou seja os produzidos junto do cátodo, 
os quais não sofreram perdas por colisão no espaço 
escuro, os quais são bastante lentos em virtude das 
colisões não elásticas que sofrerem, São estes últimos 
que produzem a zona luminosa negativa e não os do 
primeiro grupo que permitem explicar a existência dos 
máximas na zona de elevadas energias, nas respectivas 
curvas de distribuição energética. Efectivamente, quando 
a tensão aplicada aumenta, a aceleração dos electrões 
através do espaço escuro aumenta também de tal modo, 
que estes tendem cada vez mais a penetrar através do 
meio gasoso ionizado, aparecendo na região anódica. 
Este efeito de penetração é devido essencialmente ao 
facto da área de interacção com o gás ambiente dimi- 
nuir quando a energia cinética dos electrões aumenta. 
Para elevadas tensões, os electrões tendem a formar 
um feixe monoenergético, retendo esta característica 
para distâncias consideráveis através do meio gasoso. 
Se o ânodo é colocado a dada distância do cátodo, pode 
então recolher estes electrões de alta velocidade, os 
quais possuem uma energia praticamente do mesmo 
valor que a tensão de descarga, em equivalentes ev. 
De facto, observando-se por exemplo a figura 5, pode 
constatar-se que à medida que a tensão de descarga 
é aumentada, mais as curvas de distribuição de ener- 
gia se concentram em volta do valor correspondente a 
essa tensão aplicada, tornando-se o corte 
pronunciado. 


também 
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6 — IMPLICAÇÕES E CONCLUSÕES 


Os resultados encontrados relativos à distribuição 
energética de particulas negativas e nomeadamente a 
existência dos máximos na região das elevadas ener- 
gias, que confirmámos serem fundamentalmente devidos 
a electrões, têm imediatas e interessantes consequên- 
cias de grande importância, tanto na interpretação dos 
processos de nucleação de películas finas usando o 
método citado, como no conhecimento das caracteris- 
ticas das descargas luminosas em meios gasosos. Assim, 
tomando os dados publicados recentemente por vários 
autores, fixando que o uso da evaporação por impacto 
implicava na primeira fase do crescimento das películas, 
uma densidade de ilhas superior à encontrada para O 
processo de evaporação térmica no vácuo e uma maior 
facilidade de crescimento epitaxial do depósito factos 
para os quais não foi dada explicação plausível, vemos 
que vêm precisamente ao encontro dos resultados que 
conseguimos, Na realidade, a existência de feixes de 
electrões de elevada energia a bombardearem o subs- 
trato durante o crescimento da película, terá como 
primeira consequência, a criação de pontos especiais 
de nucleação, aumentando desta forma a densidade ini- 
cial de ilhas e portanto o valor da densidade de satura- 
ção de ilhas, definida como o valor a partir do qual o 
número de ilhas deixa de aumentar para dar lugar à fase 
de coalescência (ver Lewis & Campbell('). 


A elevada energia com que os electrões embatem 
no substrato, implica a alteração da superfície deste, 
provavelmente provocando através interação daquelas 
partículas com o substrato, defeitos pontuais na super- 
fície, por dissociação desta, ou obrigando os átomos 
do substrato a afastarem-se das suas posições normais 
de equilibrio na estrutura. Em qualquer dos casos, pon- 
tos especiais de nucleação serão formados, para os 
quais a energia de absorção é maior que a relativa a 
um ponto de absorção normal. O número de defeitos 
superficiais será governado pelo equilibrio estabele- 
cido entre as gamas de produção e de perda, sendo 
esta última dependente da temperatura, 


Confirmando estas hipóteses, podemos referir os 
trabalhos de Stirland e Palmberg . Stirland, idealizou um 
sistema para deposição de ouro sobre cloreto de sódio, 
usando evaporação térmica no vácuo, o qual incluia um 
substrato composto por duas partes colocadas lado a 
lado, mas de tal modo que uma delas pudesse ser bom- 
bardeada por um feixe de electrões durante qualquer 
fase da deposição. Assim, todos os factores interve- 
nientes no processo seriam iguais para as duas partes 
do substrato, excepto que uma teria a possibilidade de 
ficar submetida a bombardeamento electrónico. 


Os resultados obtidos, mostraram que a presença 
de electrões durante a formação da película, influen- 
ciava quer o processo de nucleação e crescimento desta, 
quer o seu crescimento epitaxial, mostrando ainda que 
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um depósito bombardeado, se apresente mais continuo, 
após conveniente duração da experiência, 


Hill encontrou efeitos similares, ao depositar ouro 
sobre substratos de vidro, na presença de electrões 
incidentes. Por sua vez, Palmberg através de experiên- 
cias semelhantes às de Stirland, comprovou experimen- 
talmente que a presença de electrões durante a forma- 
ção do depósito, induzia o seu crescimento epitaxial 
À este respeito, está ainda por esclarecer se o bombar- 
deamento electrónico produz o crescimento epitaxial 
através dum mecanismo que melhore a orientação dos 
núcleos, ou através dum mecanismo que aumente a 
densidade dos núcleos reponsáveis pelo crescimento 
epitaxial, cuja coalescência com núcleos não orientados, 
produza uma recristalização destes. 


Em resumo, sendo aparente que o bombardeamento 
electrónico aumenta a densidade inicial de «ilhas» ou 
núcleos estáveis, não está ainda decidido se esta con- 
dição ou outra são suficientes para explicar as diferen- 
cas no crescimento epitaxial, Acrescente-se que Layton 
encontrou efeitos similares com iões positivos a bombar- 
dearem películas finais de ouro. 


No que respeita às implicações dos nossos resul- 
tados nas características das descargas luminosas em 
meios gasosos, pretendemos apenas focar que a exis- 
tência de electrões de elevada energia junto do ânodo, 
envolve um percurso livre médio bastante grande para 
aquelas particulas, o que por sua vez obriga a que a 
área de colisão (elástica e não elásticajneste caso 
relativa ao AÁrgon tenha uma variação em função da 
energia dos electrões que atravessam aquele meio um 
pouco diferente da que nos tem sido dado apreciar 
através dos dados existentes, diminuindo muito rapida- 
mente quando a energia dos electrões aumenta, a partir 
de certo valor, mas seguindo no entanto a mesma forma 
de variação já exposta anteriormente. 


Como observação final, registe-se que no seu 
aspecto geral, os estudos teóricos relativos ao fenómeno 
da pulverização catódica estão ainda numa fase pratica- 
mente elementar, pois se a influência de alguns parã- 
metros está amplamente conhecida, o desconhecimento 
de outros não permite um completo domínio dos fenó- 
menos em jogo. 


Os resultados relativos a partículas positivas mos- 
traram-se de consequências menos interessantes e simul- 
tâneamente de mais dificil obtenção, pois as correntes 
em jogo são de valor bastante inferior às relativas a 
energias negativas. Estes resultados, farão parte duma 
futura exposição, em que se apreciará também os efeitos 
devidos a particulas neutras. Actualmente, o nosso 
principal foco de interesse, concentra-se no estudo das 
características de resistividade de películas finais 
submetidas a bombardeamento electrónico durante a 
deposição. 
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"El estano y sus aplicaciones” n.º 101 


LA TIPOGRAFIA MODERNA SOBRE ESTANO 


Una de las numerosas caracteristicas que han con- 
tribuído a que la hojalata sea un material primordial de 
embalaje es su idoneidad como superficie tipográfica, 
permitiendo la obtención de ilustraciones brillantes, con 
reproducción fiel y precisión de detalles, Además de 
su empleo en envases para alimentos, la hojalata impresa 
encuentra aplicación en los artículos denominados de 
«uso general» tales como cajas decorativas para gallatas 
y bomboneria, papeleras y bandejas así como calas her- 
méticas para baterias. El impresor de trabajos por encar- 
go ha de ocuparse de todas estas aplicaciones y en un 
artículo del último número de «El Estaho y sus Aplicacio- 
nesy se describen las actividades del principal impresor 
de hojalata de Inglaterra. Se bosquejan los fundamentos 
de la tipografia sobre hojalata y las innovaciones tales 
como los métodos electrónicos del tratamiento de tra- 
bajos artísticos y los procesos tipográficos a tres colores. 
Estos adelantos en la tecnologia tipográfica, junto con 
las propiedades intrinsecas del metal, mantiegnen en 
vanguardia a los productos de hojalata para el embalaje 
y otras aplicaciones. 


EQUIPOS DE DESTILACIÓN CON REVESTIMIENTO DE 
ESTANO 


En muchas aplicaciones, tanto médicas como indus- 
triales, existe la necesidad de agua de elevada pureza. 
Una de las formas más eficaces de conseguiria es 
mediante destilación pero los materiales de construcción 
han de ser de tal tipo que no contaminen el agua. El 
astafo ha demostrado su superioridad en coste, duración 
Y no toxicidad con respecto a otros materiales competi- 
tivos en dicha aplicación, En «El Estao y sus Aplica- 
ciones» n.º 101, un directivo de un importante fabricante 
de equipos de destilaciôn, con revestimiento de estao, 
describe la fabricación de alambiques y de otros equipos 
y el papel vital desempefado por el estaho. 


LOS PRODUCTOS QUÍMICOS DE ORGANOESTANO EN 
LA INDUSTRIA DE PLÁSTICOS 


Los compuestos de diorganoestano en los cuales el 
átomo de estaho se encuentra directamente enlazado 
con dos átomos de carbono, tienen cierto número de 
aplicaciones importantes en determinados plásticos de 
forma que su desarrollo va unido al de la industria ce 
plásticos. En un articulo del presente número se des- 
criben recientes progresos tecnológicos, tanto en los 
compuestos de diorganoestaão como en determinados 
plásticos. Los estabilizadores de organoestaão son los 
de mayor eficacia conocida para mantener la transpa- 
rencia y resistencia a los golpes del CPV durante su 
elaboración y posterior servicio. Se resenan los adelantos 
alcanzados con estos estabilizadores, incluyendo el desar- 
rollo de los tipos con «mezcla de metaly en los cuales 
los compuestos de organocestafo se combinan con sales 
de bario o calcio y los nuevos estabilizadores de dime- 
tilestafo que, según se reivindica, son estabilizadores 
térmicos de gran eficacia, Los adelantos que se describen 
en la tecnologia de los plásticos comprenden el CPV en 
aplicaciones de construcción; las espumas de poliuretano 
catalizadas con organoestafo para uso en mobiliario y 
suelos sin uniones; y la producción de zapatos con 
ampleo de moldes de silicona, que se vulcanizan a tempe- 
ratura ambiente por medio de la adicióôn de un agente 
secante de organoestafio. 


Este último número de «El Estafio y sus àAplica- 
ciones» contiegne asimismo un informe sobre el 500 
Aniversario ce la Worshipful Company of Pewterers, 
noticias del Instituto, incluyendo la inauguración oficial 
en Sydney, Nueva Gales del Sur, del Australian Tin Infor- 
mation Centre, noticias sobre próximas reuniones, criti- 
cas de libros y resumen de recientes publicaciones del 
Instituto, 
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Dispersão de poluentes na atmosfera 


Em Setembro de 1972, o Núcleo de Estudos de 
Engenharia Mecânica (|. A. C. |) e a Divisão de Termo- 
dinâmica Aplicada (|. 5. T.) organizaram no |, 38. T. um 
Curso de Extensão Universitária sobre a Dispersão de 
Poluentes na Atmosfera. À lição inaugural desse curso 
foi já publicada dados os aspectos muito gerais e pouco 
especializados nelas tratados. Segue-se agora o texta 
das lições proferidas pelo Prof. J. Mathieu. Estas lições 
formam, em si mesmo, um todo coerente e correspon- 
dem a uma contribuição original e significativa para o 
estudo da difusão na camada limite atmosférica. 

Posteriormente serão publicadas as lições sobre 
«Camada Limite Atmosféricas e «Métodos Numéricos 
de Previsão», A primeira destina-se apenas a introduzir 
os conceitos fundamentais associados à camada limite 
atmosférica nos aspectos em que diferem da habitual 
camada limite hidrodinâmica. As outras tem como fina- 
lidade uma mais ampla divulgação do trabalho de inves- 
tigação realizado no Núcleo de Estudos de Engenharia 
Mecânica relativamente aos métodos numéricos de inte- 


gração das equações fundamentais. Espera-se deste 
modo que o conjunto de publicações que agora se inicia 
B que constituem o prolongamento natural de anteriores 
trabalhos sobre «Modelos de Turbulência e sua Verifica- 
ção Experimetals também publicados pela Técinca, for- 
mem uma base sólida e cientificamente evoluída para 
a análise dos problemas postos pela. dispersão de polu- 
entes 

É numa óptica científica que desejariamos ver 
tratados os problemas de poluição que assumem entre 
nós aspectos bem mais alarmantes do que por vezes 
se supõe sem apresentarem contudo as catrastóficas 
implicações que certa imprensa insiste em nos divulgar 
a algumas actividades gostam de canalizar em beneficio 
próprio, 

Abrindo as suas páginas a trabalhos de indole 
científica sobre poluição e fomentando o interesse pelo 
tema, a TÉCNICA presta um serviço relativo que 
muito nos apraz registar e a que procuramos corres- 
ponder com as publicações que se seguem, 


Jd. J. Delgado Domingos 


C. D. U, 533.15:551.554:628.52 


Diffusion dans les couches limites stratifices 


RESUMO 


Estuda-se a estrutura das baixas camadas da atmos- 
fera usando uma atmosfera padrão em referência à qual 
é analisado o escoamento real sob a forma de uma 
perturbação. 

Estabelecem-se os perfis de velocidade e tempera- 
tura médias na camada limite atmosférica para várias 
condições de estabilidade. 

Faz-se o estudo da difusão turbulenta pela ana- 
logia molecular e pela teoria estatística sendo analisadas 
as dificulades decorrentes da não linearidade do fenó 
meno. Mostra-se como uma perturbação térmica pode 
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J. MATHIEU 


ABSTRACT 


The structure of the atmospheric boundary layer 
is studied by means of a «deviation flow» from the 
standard atmosphere. 

Mean velocity and temperature profiles are stablis- 
hed for several stability conditions. 

Turbulent diffusion is studied both through a gra- 
dient transport approach and through a statistical the- 
oryp. Difficulties due to the non linearity of the pheno- 
mena are analysed. lt is shown how a destabilizing 
thermal action can speed up turbulent diftusion. 

Measurements in wind tunnel and its limitations 
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acelerar o processo da difusão turbulenta ou pode blo- 
quea-lo consoante o tipo de estabilidade da atmostera 

Referem-se os ensaios efectuados em tunel aero 
dinâmico, as limitações e vantagens em relação às 
medidas «in situp, As dificuldades de transposição 
dos resultados do tunel para a atmosfera são devidas 
em primeiro lugar à pequenez das escalas de tempo 
e comprimento no túnel e à dificuldade em simular regi- 
mens com instabilidade significativa, 

Apresentam-se vias para melhorar estes aspectos. 

A comodidade e economia das medições feitas 
em tunel justificam a sua utilização e desenvolvimento 
mesmo que sejam aplicadas apenas para o estudo de 
tendências ou efeitos globais, 


CHAPITRE | 
Structure de la basse atmosphere 


Nous consacrerons l['ensembla de ce chapitre à 
Vétude de la structure des basses couches de l“atmos- 
phêre, Référence sera alors faite à une atmosphêre type. 
Les déviations par rapport à ce modéle idéal seront 
ensuite analysées. Pour des raisons de simple commo- 
dité cet «écoulement de déviation» sera décomposé en 
écoulement moyen et en écoulement turbulent, Afin de 
donner les équations contrôlant I'évolution de chacun 
de ces mouvements élémentaires nous établirons d'abord 
les équations instantanées du mouvement de déviation. 


1.1 Equations instantanées du mouvement de 
déviation 


Les équations du mouvement d'un fluide visqueux 
compressible soumis à [action de la pesanteur sont: 


lFéquation de quantité de mouvement: 


GER Qua E a 


E ou, u, = O 
“OL Dx Ox 


| équation d'ênergie prise ici sous la forme 
d'une équation d'évolution de l|'entropie: 


os os 

e a = 
PT(= 55) 
VU;| 4 UU; 
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and advantages are compared with «in situn measure- 
ments. Difficulties in transposig the wind tunnel measu- 
rements to full scale are mainly due to the small wind 
tunnel time and length scales and to the difficulty of 
simulating fairlp unstable flows, Several ways of impro- 
ving the measurements are presented, 

Economy and convenience of wind tunnel measu- 
rements justify its utilization and development even if 
our purpose is to seek tendencies or global effects. 


On voit qu au cours du mouvement | entropie de la 
particule fluide est accrue du fait des processus dissi- 
patifs et altérée par les éêchanges thermiques, accrois- 
sement et altération restant pondérés par le facteur 1/T. 


— l['éêquation d'état, On supposera que l'atmos- 
phêre est constitués par un gaz parfait de com- 
position constante obéissant à la loi: 


E a RT 
Ê 
Quelque soit l'état de latmosphêre nous ferons les 
hypothêses simplificatrices sulvantes;: 


a) les vitesses sont faibles devant la vitesse du 
son, 


b) les variations de pression, masse volumique el 
température demeurent relativement faibles en 
ce sens que les écarts de ces quantités par 
rapport à un état standard à définir demeurent 
petites comparées aux quantités elles-mêmes, 


c) à une hauteur donnée, les fluctuations de pres- 
sion sont sans influence sur la densité, le gaz 
adoptant la pression relative à cette altitude, 


d) les variations de à et yu avec la température 
peuvent être négligées. 


Puisque nous venons d'introduire la notion d'écart 
vis-à-vis dun modêle standart il convient tout d'abord 
d'en définir les propriétés. L'atmosphére standard est 
au repos et vérifie les relations suivantes faciles à dé- 
duire du groupe précédent dont se trouve évidement 
exclue I'équation de continuité, 


à. DP dê, OR ma 
E dx dl Spas ST 


BRT 
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f 4 
fa 
Nous écrirons encore la variation d entropie 


Quant à l'équation de conservation de la masse on 
voit en la développant qu'elle se met sous la forme 


Se 2 Pai) sl CAS sê, dá 
RT E: (7) +& + adP Cher Po) “(P +p') DEU it Meto o 

à ud ART atá ST RT AP a í Os 0%; 

e pum a et que [$)= POP P 


+ 


GU C, aT qaP aP Compte tenu de la permanence dans le temps de 
T P [Fatmosphêre standard et des variations spatiales relati- 

vement modérées de p. elle se réduit finalement à 

Dans le cas du modêle standard nous la suppose- 

rons nulla ce qui en liant les variations relatives de P à 

celles de T fixe définitivement le modéle 


fe 
Ox, 
ê a T é RaÉ 
T E, 


L'équation de quantité de mouvement prend dans 
ces conditions la forme simple: 


Dans ces conditions 


le champ thermique vertical Su. ou. A = 
sera défini par dE A =- é a -P9% 
01, - R to Bl e 3 Mais: 
03 CE Da Cp | 
49P 4 oP, 2P o 406 19P Por 
et l'on pourra écrite une relation semblable à celle Pp DR FAS 2x; Ea “Po fp st ox, 
donnéa pour la pression 
ed E. | 
a +- =R po = 2a ii BL Dans l'atmosphêre standard A Si s exprime 
É P “a en fonction de get f/p est égal « Le. on a donc: 
L'état de l'atmosphêre réel sera repéré par rapport 
à cet atmosphêre standard si bien que Ion écrira - ç , SP + 
A 
na = a E RR Sd a t 
oa O “tr q! 
e =) ) Ê : + 
Per mê J T=-1,+T | P=f+f 


L'équation de quantité de mouvement s écrit 


On supposera cependant que les écarts P, T' et p 
demeurent faibles 


OU: k TE: oU: À oP' á 
mmemmnas ET— = == —— q a sã 
OL Ox; E 'JXe Ox; ga 
Pe P . ge T a 
' Ed RM To vd 5 é forme que nous utiliserons par la suite 
Puisque les variations de masse volumique dues 


j E i En ce qui concerne l'équation d'énergie nous 
aux variations de pression sont négligées on aura ; E 

reviendrons à une formulation en T à partir de la relation 
précédement rappelée: 


t 
— f 


Il faut noter que pour conserver |'homogénéité rigoureusa dans les notations les composantes de la vitesse devralent Gtre 
notées Ut, 
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Ainsi se retrouvera mis en évidence le travail des 
forces de pression et Ion aura, compte tenu des simpli- 
fications introduites et en négligeant la chaleur dégagée 
par frottement: 


Nous écrirons alors l|'équation d'énergie pour les 
déviations T', p, P' comptées par rapport au modéle 
standard. Dans ce dernier les gradients de température 
seront supposés suffisamment faibles pour que leur 
influence sur les échanges soit négligée, il vient alors: 


& + e ue 2.4 — 
(e e da; LC “Ox daf 
ou a est la diffusivité thermique: = PC “ll 


n'est évidemment pas exclus de rencontrer loócalement 
dans latmosphêre réelle des gradients thermiques impor- 
tants entrainant des échanges considérables du fait du 
champ de déáviation. 


Eu égard aux conditions de définition de l'atmos- 
phêre standard 


[équation d'énergie pour la déviation T' ne retient qu un 
double phénomêne de diffusion convection associés et 
s écrit: 

> ST” ad 


a z. 
2%; Dx; 


Remarque: 


Il peut être fait usage d'une températurs poten- 
tiale définis de la façon suivante: 


T df | 
Ada E É. RP EE uu dê 
E O 
Ainsi 


E 


| Pg 


Pp5é, 


Dans le cas du modéls standard P,, pu TT. on a 
puisque 


dS=0 dO. o 
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Une atmosphêre ayant le profil de cette atmosphêre 
standard est dynamiquement stable, un volume fluida 
se déplaçant vers le haut ou vers le bas a toujours la 
densité locala variations de tempnérature et de pression 
se combinat de telle façon que Jt/93 O. La variation 
de T, en fonction de la hauteur est lige au fait qu un 
volume de fluide se déplacant verticalement est soumis 
à un refroidissement adiabatique dã à I'évolution de la 
pression, Pour un écoulement faiblement destabilisé 
les variations d entropie du mouvement de déviation 
seront donc liges à celles définies à partir de la tem- 
pérature potentiells 


1 DO A OT. o ST: 


e — .00— 


> — 
O) o3 a Eis 3 3 
1.2 Eguations du mouvement moyen de déviation 


La mouvement de déviation était défini à partir 
du mouvement réel dans |atmosphêre en prenant comme 
référence un modéle standard que nous avons défini et 
qui correspond à un systême dynamiquement stable. 
Nous supposerons sans discussion quil est possible 
da définir pour ce mouvement de déviation des valeurs 


moyennes 


si bien que l'on retrouve les liaisons traditionnelles 
entre grandeurs instantanées, moyennes et fluctuantes 


Uia U, + u; 
TO e To a 6 
Par” +p 


En supposant que les équations instantanées du 
mouvement de déviation soient vérifiges même lorsque 
celui-ci présente un caractêre turbulent on a en prenant 
les moyennes et compte tenu de quelques arrangements 
classiques: 


DU OU dum 4oP OU, 
dE ” OX, "Om, a Eos das 

«2T8, 
VU, » 
DL; E 
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2 Tr —g ams 
| AU DI, 000  q'T 


» “ox Vxj der 


La doubla nature des contraintes laminaire et turbu- 
lento peut être prise en compte par un tenseur unique du 
second ordre groupant éventuellement les actions norma- 
les ou de pression: 


o 3 — — o RR e 
fi o e E o A UU. O — [ef 
dx; Ê du jo e + Dx; é) ER e 


Le tenseur prenant en compte la participation du 


mouvement turbulent aux actions de contact = p, u, u; 
est gênéralement appelé tenseur de Reynolds. 


Le double caractêre des échanges thermiques donne 
ligu à un regroupement analogue si l'on introduit le 
vecteur de composants H; val que 


Nous préciserons certains aspects du problême en 
considérant un écoulement permanent dont [unique 
composante de la vitesse moyenne est U,, les proprietés 
demeurant homogênes dans le plan x, x:. On a alors: 


U: =V [x,) Bi 


Les équations de quantité de mouvement et d'êner- 
gis peuvent s'écrire de la maniêre suivante: 


Va == 
oP VV ae 
— (amç as = — Ye DER PLA 
ax dxr DX; 1h 
A La ! 
O (0, ) x oT 
DAs Dxy 
soit en faisant successivementi = 1,2,3: 


Ceia RE pl 5 fixa, v2Y 19P 


Dx; Poxa oxt DER sai “Ed, 
2 oP” 
am E qe 
di bg 4 a E E em 
Du 40,9 mem o. AA += E 
Dx, Ps Ds T “Que P. To 
) Dr — Tr! 
mg a Dom [Monde Da 
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Introduisons la double notation 


-Pma, + ra = És 


- LC ,0 4 97 =H 
Dx, 


et considérons le développement d'une couche limite 
la long d'une paroi plane. Si I'on postula [existence 
d'une zone à contrainte tangentiells et flux constant, 
on a d'une part 


Eis > fu Lp =U, - ut 
Dx, Ap, Fo 


u, désignant la vitesse de frottement, 


et dautre part: 


das) 
DL, 


= > H,=[ CU, T = ut 


T, désignant la température de frottement: 


ja E 
POL, fel 


La température de frottement rapporte comme on 
la voit le flux thermique à un flux convectif dont la 
vitesse caractéristique serait précisément la vitesse de 
frottement 


1.3 Equations du mouvement fluctuant de déviation 


Les équations du mouvement instantané peuvent 
être écrites en tenant compte pour le mouvement de 
déviation de la décomposition adoptée entre grandeurs 
moyennes et grandeurs fluctuantes, si de ces relations 
on soustrait les équations du mouvement moyen déjã 
établies on met alors en évidence |'histoire des différents 
parametres turbulents du systême. C'est ainsi que l'on a 
an désignant par «k» |'indice de sommation 


pls = 4 , 
(4) ão + Mg + Micas aC dai 
L DX O Va, “ox, 


Err 
Po Di Dão A 
Os o 
Ux, 
vo OT 0% 00 ,% do 
Peças Rc bs ft a 
OL DX 2 DX, DX, DXE 


ao 


Examinons tout d'abord le bilan d'énergie méca- 
nique. |l est commode d'opérer de façon classique en 
passant par |'intermédiaire de I'équation des contraintes. 
Nous considérons lI'équation (1) et son homologues 
obtenu en permutant les indices i et | 


(29 Eus - Pior +U fd Peão Rca 
e Ec op + Y O ui Je $ 0d 
o 0X, dxp q 


Ces équations sont alors multipliges par u et u 
respectivement puis ajoutées, on a aprês avoir une 
moyenne sur une série d'épreuves et compte tenu de 
quelquer réarrangements: 


ob dg 2x 2x, 


La corrélation u, u, d'une structure turbulante que 
l'on suppose suivre la mouvement moyen est affectéa: 


par la mouvement moven lui-smême,le terme 


praenant en charge ce couplage, 


— par le champ turbulent agissant sur son propre 
développement*, 


— par la viscosité dont le double rôle sera rappelé 
par la suite, 


— par la champ thermique associé, 
| est facile de procéder à un bilan énergétique 


portant cette fois sur le mouvement fluctuant, faisons 
dans l'expression précédente i = j et sommons i, on 


Le rôóle de chacun des termes est facile à com- 
prendra eu égard aux commentaires donnés pour I'équa- 
tion des corrélations. Toutefois, pour mieux saisir le 
rôle de la viscosité, il y aurait lieu d'appliquer cette 
relation à un certain domaine fluide, on comprendrait 
alors l'intérêt de distinguer ce qui se passe à la fron- 
tire de ce domaine et en tout point interne. Ceci revien- 
drait à mettre en évidence lopérateur divergenca dans 


Vexpression FERE SC On écrirait alors: 
Me Ge ——— “ 
pia gy é [4 é gif, À 
Da, Dar NS * Dad 
a 


ter groupe 


a 


2ême groupe 


La ter groupe de termes rend compte du transfert 
d'énergia par viscosité. Celle-ci joue un rôle tout à fait 
analogue à celui rencontré dans les équations du tour- 
billon oú le rotationnel se trouve diffusé par Vintermé- 
diaira de ». Ce movyen de «porter» en tout point du 
fluide une information qui se trouve être dans les équa- 
tions d'Helmoltz le rotationnel et qui est ici ['énergis 
cinétique turbulente q” est dans le cas d'un mouvement 
turbulent particulisrement peu efficace comparé à celui 
qui résulte de l'agitation même du fluide; par la suite, 
cette action visqueuse sera négligée. Le second groupe 
rend compte du róle dissipatif de la viscosité qui lui 


se manifeste en tout point du volume. Ce rôle est 
fondamental, une distribution plus ou moins continue 


de puits d'énergie s organise dans le domaine particu- 
liBgrement alimenté par les petites structures; on sait 
que ce taux d'énergie à dissiper est fixé par ailleurs, 
assurêment par une classe de tourbillons assez étroi- 
tement corrélés au champ moven qui assure efficace- 
ment les premiers stades de leur développement**. 


Parallâlement à ce qui a été fait pour |'ênergie 
cinétique, nous écrirons maitenant |'équation donnant 
Hevolution du carré moyen des fluctuations de tempéra- 
ure. Nous avons précisé au début du paragraphe le 
groupe des relations de base. La dernidre étant multi- 
pliée par 6, et une moyenne effectuée on trouve en 


aura en possant posant: 
du; u, 
* Mis sous la forme y,  ——— on voit mieux [aspect convectif de ce terme la mouvement turbulent se váhiculant 
t) Xj p 
lui-même. 
** Ceci laisserait à penser que les axes des trourbillons qui prétêvent l'énergie de I'écoulement moven sont par exemple 


dirigés suivant laxe principal de la quadrigue de déformation qui assure au mieux leur croissance, l'axe de dilatation par 


conséquent. 
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[e 
= 

= + D= 5 DT 2 — o 

a l a dm k Dx, Da, “ox 


Cetta équation est à comparer à celle obtenue 


pour q. L'énergie thermique en suivant le mouvement 
moyen est affectéc: 


—- par l'inhomogénéité spaciale du champ de 
déviation; le terme que l'on peut interpréáter 
comme une production 


prend en charge le couplage entre champ ther- 
mique moyen et champ turbulent, 


— par le champ turbulent agissant sur son propre 
développement: 


— par la conductibilité thermique. Le terme 1028. 
dxt 
k 


qui la prend en charge est souvent supposé 
peu différent du terme essentiellement négátif 


E NO Se) ce qui revient à négliger le 
A Pl q 28 c'est-á-dire à suppo- 
terme DXn ox, 


ser la turbulence thermique homogêne du point de vue 
de la conduction. 


Si 


he (a eo To + 50) =0 


on a bien en effet 


Remarques 


1.) Dans cette équation la pression ne figure 
d aucune maniêre ce qui se conçoit si Ion tient compta 


de son rôle répartiteur que nous préciserons par la 
suite et du caractêre scalaire de la grandeur traitée” 


dx. 


2.) Nous avons interprété le terme 


comme une production d énergie turbulente. Pour 
sen assurer il conviendrait de procéder comme on 
a coutume de le faire lorsque l'on cherche à donner 


un sens au teme aqu, Pis,  c'estâ-dire en 


examinant |'équation d'évolution des températures moyen- 
nes du champ de déviation: 


EMO au É. é eu 4 9H 
E x, xe dl DX Cp DX 


avec toujours 


F 


or PRA 
Ho=- AS +C Ou 
Xu T' 


En effet multiplions cette équation par T', elle 


devient:; 
DPM ar" dart hM 
SE 3 ERA at 2 Co Ox 


dont l'action s'exerce en tout point mais il est affecté 
d'un signe plus autrement dit s'il constitue un apport 
pour le champ fluctuant thermique & il constituera un 
prélêvement vis-á-vis du champ des températures 


moyennes T', 


Les équations des corrélations vitesse température 
s'écrivent en procédant de la même maniêre. On recon- 
nait dans les deux équations suivantes celles qui nous 
ont servi à former |'équation d'évolution des corrélations 


uv u; et 6º multipliées respectivement par u et1 puis 


movyennées., 


*Une comparaison avec l'équation du tourbillon est souvent de ce fait tentée. 
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E dj “ k OXy 
Rn o ºP ++ q du. A. “SE 
pd; o E 
Te + dl piA E ST +U, PP aa je LAS Mico 
E Xp k JL, 
E: 
ro E A 


Opérant à la fois une sommation et un regrou- 
pement des termes on trouve: 


D Bis o Me CÊM gm QUI e OT 

a = á 

= qua! 00 « 8:68 » gm Ve 
ox PCA att o * "OS tia 


L'évolution des corrélations vitesse-température en 
suivant le mouvement moyen est affectée: 


— par l'inhomogénéité spatiale du champ moyen 
cinématique dU, / dx, et thermique dT'/ dx, cia 
prorata des corrélations déjã existantes qu, et 


u; u, respectivement, 


-—- nar le champ turbulent agissant sur son propre 
développement, 


— par les corrélations pression-température; si la 
pression disparait de l'équation fixant l'évolu- 
tion des fluctuations de température du fail! 
du caractêre êéminemment scalaire de cette quan- 
tité, elles s'introduisent au niveau des corréla 
tions température-vitesse, 


| 9! Bo 


— par laction des forces d Archimêde —, 


—— par le rôle conjugué de la viscosité et de la con- 
ductibilité; le nombre de Prandtl jouant le róle 
de coefficient de pondération 


. , au ha , | ! 
W os da -v jo E a 


Des == 
Dep Dx oxt SG? ok] 


Ces termes,dans lesquels figurent des dérivêes 
secondes sont évidemment actifs sur les fines strue- 
tures du champ. Viscosité d'une part, conduectibilité 
thermique d'autre part tendent à atténuer les variations 
brusques de vitesse et de température, cet effacement 
partial n étant pas nécessairement aussi actif dans les 
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deux cas, au moins lorsque le nombre de Prandtl est 
différent de 1, SUNYACH et MATHIEU [22] puls 
CORRSIN et KUO [33] signalent dans l'air que les 
structures turbulentes thermiques paraissent réparties de 
tacon plus homogêne que les structures turbulentes 
cinématiques. 


Examinons & nouveau le cas d'un écoulement 
permanent dirigé sulvant x, et possédant des propriétés 
d'homogénéitée dans le plan x, x. On avait posé 


E, = Víx,) dis 


(c'est des modéles simples de ce genre auxquels nous 
nous référerons ultérigurement) et écrit les équations 
du mouvement movyen de déviation. Áfin de compléter 
nos connaissances sur le problême d'écoulements quasi 
paralléles, revenons aux équations des corrélations dou- 
bles en un point dans lesquelles nous ferons i = j. 
Dans un premier temps, nous ne procêderons à aucune 
sommation sur cet indice répété de façon à saisir com- 
ment est alimenté le mouvement turbulent du fait du 
double processus mécanique et thermique. On a succes- 
sivement et en négligeant la diffusion lige à v: 


Es — > as 2 
det DU, 8 00 apa 
o: s DX [88 dx, dE) 

Õ= o AA, LA er op - EV du, 
Ox, & *Dk, DX, 


On y retrouve une situation analogue à celle 
rencontrée dans une couche limite non chauftée. Négli- 
geant l'action thermique on voit en effet que seule la 
composante u, est alimentée par [écoulement moyen. 
C'est alors le champ de pression qui en tout point 
redistribue entre les différentes composantes | énergie 
fournia à la seule composante u,. Dans le cas qui nous 
occupe deux comnosantes sont alimentées par le champ 
movyen, non seulement la composante u, bénéficie d'un 
apport du fait du gradient de vitesse, mais aussi la 
composante us par le truchement du champ thermique 
associé. L'action de ce dernier, sur les termes prenant en 
charge la diffusion, se décelle à travers u, et Ion retrou- 
vera en effet, aprês sommation dans l'équation d 'éênergie 


un terme u. ql. Nous vérifierons par la suite que 
Vinfluence du champ thermique sur la diffusion est três 
afficace, 

Pour fixer les idées, il peut être intéressant de 
préciser la situation rencontrée dans le cas d'un écou- 
lement destabilisé par chauffage d'une paroi inférieure, 


on a alors les profils cinêmatiques et thermiques suivants. 
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profil thermio Um 


profil cinômatigque 


PPFFFFFFFFFERFI 


Du point de vue moyenne statistique à des u; posi- 


tifs sont associéês des u, négalifs de sorte que 


at 


La même raisonnement montre quavec le profil 


thermique ci-contre tu, > O de sorte que les deux 
tarmes producteurs dans les équations précédentes con- 
juguent leur action. Si la paroi inférieure est refroidie 
ves Jeux mêmes termes tendent alors à se neutraliser 

Procêdons  maintenant à une sommation sur 
Cindice |: 


Va e Pp ei du 
o mim ND fe(ato E) o(s Dans 


Nous poserons par la suita de façon à allêger 
Vécriture: 


Avec les mêmes hypothêses [équatior qui contróla 


les fluctuations thermiques 9-2 s'éerit: 
Te E RP “aum dot 
O = — | Bu — (4) +a[5) 
a! 3) * Dx, (943) dx, 
vo 
avec e 
o (35) 


Dans le cas qui nous occupe c'est la corrélation 
9 us qui présente le plus d'intérât, elle s écrit: 


mais il est également possible de faire apparaitre |his- 
toire des corrélations tu, et I'on trouve alors: 
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On voit intervenir sur chacune de ces corrélations 
qui contrôle respectivement les échanges thermiques 
transversaux et longitudinaux les gradients de tempéra- 
ture et de vitesse moyennes tandis que laction d'Archi- 
mêde ne se manifeste qu'au niveau de I'équation réglant 
les flux transversaux. 


L'examen de ces équations suggêre immédiatement 
da comparer la production thermique à la cinématique. 
Ainsi introdui-on le nombre de Richardson 


3 u, 6 

— i b 

E + 
a ox, 


R 


Si l'on revient à l'exemple précédement traité on 
constate que le régime instable pour lequel les forces 
de poussés d'Archiméde tendent à accroitre l'énergie 
turbulente correspond à des nombres de Richardson 
négatifs (0 > O; u u < O; dV/dx, > O). Le 
régime stable correspond bien évidement à des rombres 
de Richardson positiís. 


Dans le cas ou I'on admet un liaison simple entre 
corrélations et champ moyen on a en posant 


souvent les météorologistes supposan' E. Er c'est-à- 
-dire I'égalité des deux codfficients de uitusion ciné- 
matique et thermique préférent ne considérer que la 
forme simplifige du nombre de Richardson 


plus facile à atteindre expêrimentalement et appelée 
nombre de Richardson de gradient. 


L'examen des mêmes équations nous conduit 
d'emblée à introduire une comparaison entre production 
et dissipation en posant 
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Ces arrangemens formels coduisent au bilan d'êner- 
gétique suivant: 


la dernier terme prenant en charge, au point considere, 
l'énergie importée ou exportée, L'aspect énergétique du 
problême revient donc a examiner le signe de l'expression 


(1-8) 


ou D est un nombre essentiellement positif et R un 
nombre négatif (cas instable) ou positif (cas stable). On 
conçoit alors comment le phénomêne thermique peut 
alimenter ou modêrer le mouvement turbulent, Une 
analyse des conditions ou stabilité faite par Taylor 
(1931) montre que pour R, > 0,25 aucune turbulence 
n'existe en l'absence de perturbations extérieures impor- 
tantes, En soufferie, des mesures de turbulence ont été 
faites par ELLISON (1962) jusqu'ã des valeurs deR, 
voisines de 1 encore qu'il soit difficile d'affirmer que 
ladite turbulence ne soit pas en déclin. Dans l“atmos- 
phêre par contre il semble bien qu'aucune turbulence 
n'ait pu être déceléa dês que R, excêde la valeur 0,2. 


1.4 Régimes d'écoulement 
stratifiées 


des couches limites 


Les indications précédentes permettent donc un 
classement des couches stratifiges et certains rappro- 
chements avec les modéles homogênes étudiés. 


Conditions stables: 


En dehors du cas exceptionnel ou l'atmosphêre 
serait complêtement calme on trouve au voisinage du sol 
une région ou R. est faiblement positif. L'effet du champ 
thermique n'est cependant pas suffisament intense pour 
n«étouffer» tout mouvement d'agitation. Dans cette zone 
ou la tersion demeure sensiblement constante et égale 
à la tension pariétale on peut considérer que comme 
dans une couche limite ordinaire production et dissipa- 
tion contrôlent à peu prês seules les mécanismes tur- 
bulents. 


Lorsque l'on gagne les couches plus élevéesR, croit 
et la turbulence cinématique est en que!que sorte absor- 
bée. Cette évolution de la stabilité avec l'altitude est 
rendue possible du fait que |'énergis cinétique turbulente 
a une double origine: une origine cinématique pour la 
quelle gradient de vitesse moyenne et contraintes sont 
associés, et une origine thermique ne mettant en cause 
aucun gradiant moyen,. 


Conditions instables: 


Nous plaçant dans le cadre des approximations 
météorologiques classiques nous tenterons de distinguer 
4. régions: 


| —Du fait de I'évolution relative des 2 productions 
turbulentes (cinématique et thermique) on doit s atten- 
dre à rencontrer au voisionge de toute paroi lã ou les 
gradients de vitesse sont intenses une région x. < 10m 
ou domine l'aspect cinématique. 


Les caractéristiques de I'écoulement sont alors cel- 
les relevées sur des couches limites de soufflerie, Cette 
région de prévalance est encore appellêe région de 
convection forcée. Les deux phéromênes cinématiques 
et thermiques s'y trouvent totalement découplés, le 
premier fixant les modalités de fonctionnement du 
second. Le nombre de Richardson est tel que | R, | <]0,03. 


| — Au delã de cette zone existe une région située 
dans la gamme des altitudes suivantes 10 < x, < 100 m 
< IR; < 0,1, pour laquelle 

le transfert vertical de quantité de mouvement et de 
chaleur est essentielement dô à la convection thermique. 
Ne retenant que cet aspect du problême/ elle est le 
plus souvent désignée par«région de convection libren**, 
appellation quelque peu abusive si l'on se rend compte 
que la convection libre proprement dite n'apparit que 
lorsque le nombre de Richardson est égal à 1*. Plus 
précisément, disons que les effets thermiques y créent 
des mouvements verticaux à grande échelle bien corréllés 
avec la température dont le champ est déterminant quant 
à leur formation. Ces structures turbulentes particúliêre- 
ment efficaces sur le plan des échanges ont des dimen- 
sions comparables à celle de la zone instable, du point 
de vue spectral elles correspondent à la partie basse fré- 


* Pour cette valeur énergie mécanique et thermique se produisent au même taux. 


** Nous préférerions l'appellation région de convection mixte. 
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quence du spectre se différentiant ainsi assez nettement 
des structures d'origine mécanique à fréquences moyen- 
nes et élevées. Notons encore que |'anisotropie de cette 
agitation d'origine thermique est três marquée, en tous 
cas davantage que celle d'origine mécanique. Bref cette 
région est caractérisée par la présence de larges struc- 
tures convectives d'origine thermique assurant |'essentiel 
des échanges (mécaniques et thermiques). Cependant 
si de ce point de vue on assiste à une prédominance 
de l'aspect thermique, cela ne signifie pas pour autant 
que la quantité d'énergie mécanique soit três faible 
vis-á-vis de celle li£e à ces mouvements convectifs de 
vaste amplitude, Compte tenu de ce quasi-découplage, 
Lumley et Panofsky acceptent sur le plan thermique 
lidés d'un bilan turbulent donné par les équations 
suivantes: 


ox, h, 
x, “8 + 9x, ? 4 AVEC É a 2 


La I'bre équation lie par l'intermédiaire de u, les 
mouvements convectifs verticaux à leur terme producteur 


en l'occurence la poussée d'Archimêde g/T. fu... Cette 
même composante u, intervient alors dans PSguntion 
qui fixe le niveau h des fluctuations de température. 
Notons qu'eu égard aux échelles turbulentes mises en 
jeu le fait de négliger les effets moléculaires parait 
justifié. Moins bien contrólé est l'abandon des corréla- 
tions pression vitesse dont le rôóle sur la «diffusion» 
n'a pas été non plus retenu. 


Ce modêle thermique relatif à une turbulence homo- 
gêne sera complété par un modêle mécanique lui aussi 
homogêne sans autre lien avec le précédent qu'un pro- 
cessus accéléré de diffusion turbulente permettant en 
partie à I'énergie mécanique d'émigrer vers des régions 
de la couche atmosphérique ou la dissipation se trouve 
accéléréêe par l'intermédiaire d'un processus thermique. 


ll — Plus haut entre 100 et 1000 m, se situe une 
région de transit, Le champ thermique y crée un méca- 
nisme de diffusion toujours aussi efficace par l'inter- 
médiaire des termes en U; q' mais cette fois leur relative 
homogêéneéité spatiale conduit à supposer leur divergence 
nulle. Une certaine énergie mécanique est créée du fait 
des fluctuations thermiques associêes au champ gra- 
vifiquo (IR,|> 0,1) mais un certain emiettement des 
structures turbulentes permet grosso-modo de la dissiper 
sur place. Si bien que le modêle proposé répond sur 
la plan thermique à la relation: 


* Dês que |'on veut prendre en charge des domaines importants de flulde situés en altitude, 


Iv — Au dessus de 1000m se trove une région 
d'atmosphêre stable, Les fluctuations cinématiques trans- 
portée depuis les basses couches y sont réduites par 
le jeu de poussées d'Archimêde adverses tandis que 
la viscosité affirme son action en profondeur sur les 
plus petites structures turbulentes. 


L'agitation thermique importée depuis les basses 
couches de l'atmosphêre et cella née du processus de 
réduction déjã évoquée se voient réduites finalement 
par laction de la conductibilitá thermique. Pour cette 
région, un modéla homogêne répondant aux équations 
suivantes a été donné par les auteurs déja cités: 


Õ és ua E SR + ER [5] 


dx, 7 a 
Da DT? a 
ae as Ea Ta E, 


1.5 Examen des conditions de similitude 


Revenons aux équations suivantes régissant le 


mouvement moven de déviation: 


Us Pri ge 19P' PU dT' Gt) 
dt DX, ox fot “a E” 
Cu, , 25. 45T 

jo DX. Dj “OX; 


défini par rapport à une atmosphêre neutre dont les 
propriétés sont décrites par: 


j DE g 

oiee. q! o e 
fo VAL, dá a 

4 ê ó,,=0 
Ls Cp 


les forces de Coriolis doivent 


être prises en compte, Flles apparaissent alors dans le ler membre de |'équation sous la forme Ecs Ee 
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X, 7 4000%W Dome Stable » Loma d'absaltio 


RR 
TT 


) 400< X, < 4000 mM. IR:l>o,4 


É RJ tone de ron 
100 m. T 7 T 


Zone de Pal sã Lone de comvecitom rs Te 


JO LX, 4400 MW. 0,032. |IR,|S 0,4 
Producttoy — A4O0m, 
sao Zone de con vecito) forces X, €40m. LRE|X e. 
Adoptons las grandeurs de référance suivantes: On voit apparaitre les 3 paramétres adimensionnels 
classique: 
— une vitesse moyenne U* 8 
— lg nombre de Reynolds is = qui 
— yne dimension linéaire L * 
représente le rapport des forces d'inertie aux 
— yne masse volumique P* forces de viscosité, 
— une déviation de température T'* — le nobre de Prandtl GP - LS vans leque! 
La cohérence du systême oblige à choisir pour seules se trouvent associées les propriétés 
échelle de temps physiques du fluide, 
tb qu” 
ue - le nombre de Froude t GLBT" 
a ER d/ 
et pour échelle de pression P =p U” ou B sa 1 qui compare ces même forces 


d'inertie aux forces de poussés d Archiméde., 
Opérons alors le changement de variables: 


PF. . Les météorologistes préfêrent souvent caractériser 
ao o gt Us EE EM e ; À 
diz-— 3 Str ; a, Ea la stabilité thermique de I'écoulement par la longueur 
h Es RE o dite longueur de Monin Obukhov construite sur les 3 
A? pur Tt vê ex Ag paramêtres suivants 
o q o * E “wu 'T 
f P 
il vient 
pr EEE? x 
OU! tDUL Dutus 4 DPO ADU; A TS, 9 é H, 
or "tom aa, Romi E Co o, de 
ie entEr E ” PC! É 
OT v*oT p QU 9º A T (+) 
ape & . Ei . 
ot Ox Dat RP OX; 


* Si l'on avalt conservê les forces d'inertie de Coriolis dans les équatons précédentes de façon à mieux cerner le pro- 
blême métgorologique on aurait fait apparaitre le nombre de Rossby R, = U*r/L*fo* ou q)* désigne un vecteur unitaire porté par 
la ligne des pôles. C'est le rapport entre les forces d'inertia de Coriolis et celles résultant de la non permanence et de la non 
uniformitê du champ de vitesse qui est pris en charge par R,. Pour fixer les idées disons que I'on n'aura pas à se soucier de ce 
paramétre lorsque la dimension caractéristique du site est inféricure à 15 km. 
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ou U, désigne toujours la vitesse de frottement et H, 
le flux thermique pariétal. On posse alors: 


DO o qdo 


2 Ho 
+ 4 


ka. désignant la constante de Karman. 


Introduisons la température de frottement T, talle 
que 


Hs Sc UT 


Ho vient: PoE Elis 
ed Te 


Par tradition la longueur de Monin est choisie de 
facon à être positive lorsque le flux thermique pariétal 
est négatif c'est-à-dire lorsque la couche stratifiée est 
stable, elle est négative dans le cas contraire, 


La température de référence T, peut être prise 
égale à la température moyenne dans la couche. 


Bornons-nous pour l'instant à ne considérer que 
les similitudes de Reynolds et de Froude, 


La premiére constitue le principal obstacla aux 
études sur modéle, les échelles géométriques couram- 
ment utilisées (de Iordre du 1/500ême) conduisent 
en effet à des nombres de Reynolds 10 fois plus faibles 
que sur le prototype. Cependant dans la plupart des cas 
la vent souffle sur des terrains comportant suffisamment 
d'accidents géométriques à arêtes vives pour qu'un 
régime turbulent rugueux s éttablisse. L'action de la 
viscosité sur les profils moyens de vitesse s efface et 
son rôle sur la grandeur et la maquette doit être discuté 
du seul point de vue de la turbulence: la similitude au 
niveau des petites structures dont le róle est essen- 
tiegllement dissipatif risque seule d'être gravement mise 
en défaut. 


Le nombre de Froude qui compare les forces d'iner- 
tie aux forces de poussée d'Archiméde lesquelles résul- 
tent de l'écart de température entre l'atmosphêre réelle 
et latmosphêre neutre peut facilement être la même 
dans la grandeur et dans la maquette. On réalisera 
cette égalité soit par chauffage soit par refroidissement 
du plancher choisissant par ailleurs des vitesse d'écou- 
lement appropriées généralement três faibles; dans ce 
choix la limite inférieure devra tenir compte des exi- 
gences imposées par les problêmes métrologiques. 


Nous verrons par la suite comment la longueur de 
Monin Obukhov, qui ne prend en charge que des pro- 
priétés moyennes du champ d'êcoulement s articule 
avec le nombre de Richardson, 
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1.6 Profils de vitesse et de température 


Dans tout écoulement turbulent faiblement desta- 
bilisé par chauffage les caractéristiques moyennes à la 
hauteur ,. = z dépendent comme on peut le voir d'aprês 
les êquations établies de 


do Po nã 4, eis poeta 
Le. 


Ces sept quantités étant liées à quatre dimensions 
fondamentales: lorgeur, masse, temps et température, 
nous pouvons introduire en plus ces nombres de Rey 
nolds et de Prandt| un groupement adimensional tel ir 
nombre de Froude ou mieux la rapport z/L. Dans le 
cas ou la peroi à travers laquelle se produit I'échange 
visqueuse, il faut encore ajouter aux paramêtres ci-dessus 
énoncés une hauteur z. caractéristique des rugosités 
ce qui entraine évidemment l'apparition d'un nouveau 
nombre adimensionel. 


Nous allons examiner dans ce qui suit quelles sont 
los formes de profil de vitesse et de température aux- 
quelles on doit attendre dans ce type d'écoulement. 


A — Profils de vitesse et de température en atmos- 
phêre pratiquement neutre 


L'examen de ce problêms est suffisament classique 
pour que nous ne procêdions qu 'à un bref rappel portant 
sssentiellement notre attention sur la zona turbulente 
pariêtale qui s'étend sur environ 20% de |'êpaisseur 
de la couche limite et à l'intérigur ce laquelle les varia- 
tions de la contrainte tangentielle restent faibles. 


La situation que nous allons décrire correspond 
aussi bien à celle de lamosphéêre neutre qu'au cas d'un 
ecoulement suffissamment peu chaufté pour que le phê- 
noméêre thermique n'influe pas sur le mécanique. Pour le 
problême traité il s'agit donc d'une sorte de cas limite. 


Pour faire bref 
| rapport LL — 
du, 


remarquons simplement que le 
a les dimensions d'une longueisr 
g 

Dars la couche limite pariétale qui est supposée 
ignorer la struture d'ensemble de l'écoulement et en 
particulier l'épaisseur à de la zone perturbêe la seule 
dimension à laquells ce rapport puisse etre comparé est 
la hauteur z. On a alors: 


U, | a u 
RS, RA = — = 4, —+ 


a i 

ou f, est une nouvelle constante qui peut dépendre 
soit ve la viscosité soit de la hauteur caractéristiquede 
des rugosiés. 
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Dans le cas d'une paroi lisse on retrouvera la «lo! ou B, est une nouvelle constante qui peut dépendre soil 
de paroi» oú figure la constante k de Karman de la viscosité et de la conductibilité thermique, soit 
de la hauteur caractéristique Z.. 


“ Le parallêle peut encore se poursuivre sur le plan 


eli 
E 
são 


de la production, La turbulence thermique 7 = 6/2 est 
alimentée au taux de bw dT'/dz. Si l'on suppose le flux 


Dans le cas d'une paroi rugueuse on aura, €!' 
sensiblement constant et égal à sa valeur à la paroi H, 


supposant négligeable l'action de la viscosité sur | 


profil moyen: on a alors 
u 4 ds tds (*] Wi = : Eri 
U kk dd 


L ajeiica dé Estonia WO stróu Ta Puisqu'il s'agit d'une atmosphêére neutre, la no ton 

A RES SBE  RRH MEO MHRNIDO aid cia ERES de longueur de Monin Obukhov n'a pas grand Jens 
ment intense dans la zone pariétale s exprime aisément : 

physique. Cependant on peut donner avec Monin une 


en fonction de U et z. Négligeant les variations de 7 | Sa | e 
; interprétation qui s'avére à l'usage moins fragile qu'il 
qui sera pris en tout point égal à Tp ona: 


ny parait. 
Quant à la participation du champ thermique à la 


3 Ei 

W - Ep ; av e Us, turbulence cinématique elle est égale à g/T. wi si bien 
à E nt Ep E q 

p da kg que se plaçant toujours dans la couche à tension et 


flux constant on peut écrire: 


Reprenant le même raisonnement pour le champ de 


température, il suffit de remarquer que le rapport TJdT' /d; nã , H, ui 
a lui aussi les dimensions d'une longueur. La seule Q : As lho o *T Le, MM EE 
disponible dans la couche pariétale éiant z on voit que = aa UA, 
u 
“da ha 3 
T(y)= B Te losa + B, 
Ro E ê — = ô. 
1 + E 


* Pour [ensemble de la couche limite partant ce 


| Us=tlx au, »*Jla présentation suivante peutétre retenue. Pour la zone 
turbulente pariétale on adopte la forme suivante: 


CABE) «E[B 


et pour la zone turbulente centrale une forme situant ['écoulement par rapport à lécoulement sain de vitesse Us. 


El -c(8) + 6. (218) 


à : ] [o TE 
désigne [| épaisseur de la couche. Remarquat alors que UA. et “4 Us dépendent fortement des fonctions indicées «la et fai- 


blement des fonctions indicées «Zn ceci donne au problême :m tour analogue à celui rencontré dans les méthodes de perturbation (Van 
Dycke [37]). Une tentative de raccordement des deux representations dans une bande étroite à cheval sur la limite externe de la zene 
turbulente pariétale et la limite interne de la zone turbulente centrale conduit en premiére approximation à une loi logarithmique 


SP ; au 
Au, = A +B to 5 as 
Une formulation analogue peut être donnée pour les températures 
La = C+D by 34 
F à w 


Remarquons que dans sa présentation Millikan néglige d'emblées F, et G. et son raisonnement se réduit à reconnaitre que la 
'onciton logarithmique est la seule qui permette par soustracuon de aire disparaitre des constantes. Un pourra lire également la 
ontnbution à [erude du transfery thermique de RONAT [28], 
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3 — Profils de vitesse et de température en atmos- 
phêre faiblement instable [L fini) 


Nous pourrions chercher à éclairer davantage le 
mécanisme qui nous occupe en nous appuyant sur un 
second cas limite celui de la convection libre, || s'agit 
bien en effet pour nous de définir ce que peut être une 
interaction au ler ordre entre champ thermique et ciné- 
matique. 


On peut facilement reconnaitre pour le problêmes de 
turbulence libre un processus convectif moteur liéê au 
Brme 9 Hj (PC, T.. Avec le rapport z/L ce mécanisme 
générateur (dans la forme ou il apparait pour le pro- 
blême qui nous occupe) est comparé à la production 
d'énergie turbulente à altitude z puisque 


$ u, 
os + E EE 
3) ha se 


En atmosphêre neutre une telle comparaison perd 
son sens et Ion comprend qu'alors | L | = «, 


Avec Monin et Obukhov nous introduisons alors 
la quantité 


8. ta di 
U, als 


qui pour une atmosphêre neutre doit être égale à I'unité, 
Le coefficient de viscosité turbulente cinématique rç 5'Y 
trouve liê par la relation 


RU 


Es : DUISQue: Es 


ó 


De façon analogue nous poserons 


Le coefficient de viscosité turbulente thermique 
ey se trouve lié à Up et T, par 


Dans ces conditiins on a 


R - hu,3 
. 
et ainsi 


Ss. R £r 


E 


Nous reviendrons par la suite sur la signification 
de ces deux coefficients »y et +,.; en fait ils ne se jus- 
tifientique par les commodités ou'ils introduisent dans 
bon nombre de raisonnement et plus particuliBgrement 
dans le domaine météorologique. 


C'est par rapport à ce cas idéal de l'atmosphêre 
neutre que le cas de la couche destabilisêe constitue 
une déviation du ler ordre. 


Un examen de quantités devant intervenir dans ce 
problême z, L, U, U, incite à former un gradient adimen- 
sionel de vitesse fonction de la distance réduite z/L 


df Uz ) 
DSO 2d (2) 


Multipliant ces deux quantités par kz/L il vient 


er de facon analogue 


| 
| 


R=& [%) 


Puisquil sagit de «faible déviationh on pourra 
toujours poser: 


RPa)=10E 


ou 5 désigne une constante. Alors de 


relation valable pour z >> zs (pour des raisons de 
commodité z, a été introduit comme constante d'inté- 
gration). 


Nous donnerons une variante de cette formule du 
fait que géênéralement L est inconnue les flux de chaleur 
demeurant três difficile à mesurer Dans ces conditions 
on évalua flux et contraintes à partir des gradients de 


* Seules sont prises en compte des altitudes 7 >>> kz, de telle sorte que z, ne figure pas dans ce recensement, 
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vempérature et de vitesse en introduisant les coefficients 
rg et ey, La grandeur Lº correspondante est 


E) 


2 [O Ag)T 
tea (das) 


et I'on voit qu'alors L et L' sont liés; 
JO) Er 
Ê = PR Es 


- F & Er 
Si Ion considêre EA ion plus comme une cons 
tante mais comme une fovction du rapport de z'L, alors' 


se4 (3) 


et on a en premiêre approximation 


AVBT 


La valeur de 5' devrait être constante mais Lumley 
et Panofsky signalent des divergences importantes selon 
les auteurs consultés, Elle est estimée à 4 par Deacon 
(1962) à 4,5 par Panofsky, Blackadar et Mc Vehil 
(1960) à 6 par Taylor... 
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On aura pour les profils de température une rela- 
ton analogue à la précédente en supposant ty feç peu 
différent de 1 


Rappelons que le nombre de Richardson de flux 
définl par: 


Ro. 9 Rec 
Tt (Es) 


sexprime aisément en fonction de 5º et Z/L' dans la 
zone de profils logarithmiques 


A. dp» E 
+ f3 Brg: 


Cette brêve revue nous permet donc de connaitre 
à quellas déviation des profils moyens de vitesse et de 
température |'on est en droit de sattendre lorsqu' un 
êcoulement est faiblement perturbê du fait d'échanges 
thermiques àtravers une paroi horizontale. 
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Resumo dos Artigos publicados na «Técnica» n.º 424 
Ano L— Março 1975 


C. D. U. 621.3:517.512.2 


LUIS HENRIQUE MARTINS BORGES DE ALMEIDA 
a 
GUILHERME DINIZ MORENO DA SILVA ARROZ 


ESTUDO EXPERIMENTAL DA APLICAÇÃO DO CEPSTRO (1) 
A DETECÇÃO DE ECOS E PERIODICIDADES EM SINAIS 


Técnica N.º 424 — L -—- 3-1975, p. 167 a 184 


Este artigo, que constituiu relatório de um estágio, realizado 
no Centro de Estudos de Electrônica, passa rapidamente 
em revista a técnica do Cepstro, pondo em relevo a sua 
utilidade na detecção de ecos e periodicidades e na 
demarcação, com precisão, da época da sua ocor- 
rência. Além de uma breve análise teórica da técnica 
utilizada, apresentam-se algumas simulações em computador 
a também alguns casos reais. 


C. D. U. 624.072.23 
J. 5. BRAZÃO FARINHA 
VIGAS  CONTINUAS PARCIALMENTE  ENCASTRADAS 
Técnica N.º 424 — | — 3-1975, p. 185 a 212 


Apresenta-se o formulário geral correspondente a vigas con- 
tinuas simétricas de 2 a 4 tramos, com apois extremos 
encastrados. Elaboram-se tabelas para os casos de iguais 
vãos e actuação de cargas de distribuição trapezoidal, e 
de desiguais vãos e actuação de cargas uniformemente 
distribuidas. 

Toda a dedução se baseia na aplicação do método de Cross, 
notavelmente expedito para este tipo de estruturas. 
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C. D. U. 539.232.543.42 
LEOPOLDO J. M. GUIMARÃES 


ESPECTROS DE MASSA E ENERGIA NA FORMAÇÃO DE 
PELICULAS FINAS USANDO O PROCESSO 
DE «EVAPORAÇÃO POR IMPACTO» 


Técnica N.º 424 — | — 3.1975, p. 213 a 220 


O conhecimento dos espectros de massa e de energia das 
partículas incidentes no substrato, quando da formação de 
Películas Finas é de extrema importância para a com- 
preensão dos seus processos de formação e crescimento, 
intimamente ligados às proprisgdades macroscópias que 
apresentam. Para a determinação destes espectros, cons- 
truiu-se um analizador de massa e de enegia cuja montagem 
se adaptou ao estudo do processo de «Evaporação por 
Impactos. 


C. D. U, 533 15:651.554:628.52 
J. MATHIEU 
DIFUSÃO NAS CAMADAS LIMITES ESTRATIFICADAS 
Técnica N.º 424 — L— 3-1975, p. 221 a 236 


Estuda-se a estrutura das baixas camadas da atmosfera 
usando uma atmosfera padrão em referência à qual é anali- 
sado o escoamento real sob a forma de uma perturbação. 
Estabelecem-se os perfis de velocidade e temperatura médias 
na camada limite atmosférica para várias condições de 
estabilidade. 

Faz-se o estudo da difusão turbulenta pela analogia mole- 
cular e pela teoria estatística sendo analisadas as dificul- 
dades decorrentes da mão linearidade do fenómeno. Mostra- 
-58 COMO uma perturbação térmica pode acelerar o processo 
da difusão turbulenta ou pode bloqueá-lo consoante o tipo 
de estabilidade da atmosfera. 

Referem-se os ensaios efectuados em túnel aerodinâmico, as 
limitações e vantagens em relação às medidas «in situs. 
As dificuldades de transposição dos resultados do túnel 
para a atmosfera são devidas em primeiro lugar à pequenez 
das escalas de tempo e comprimento no túnel e à dificul- 
dade em simular regimens com instabilidade significativa. 
Apresentam-se vias para melhorar estes aspectos. 

À comodidade e economia das medições feitas em túnel 
justificam a sua utilização e desenvolvimento mesmo que 
sejam aplicadas apenas para o estudo de tendências ou 
efeitos globais. 


Synopsis of articles published in «Técnica» n.º 424 


L —-March 1975 


U. D. €. 539.232:543.42 
LEOPOLDO J. M. GUIMARÃES 


ENERGY ANALYSIS OF PARTICLES ON THE SUBSTRACTE 
USING A PROCESS «CATHODIC SPUTTERING» 


Técnica N.º 424 — L — 3-1975, p. 213 a 220 


The nature of the particles incident on a substracte is od 
basic interest for on understandig of the nucleation process 
of thin films obtained by deposition. 

A sintable device has been constructed to study the nature 
and energy of those particles, using a «cathodic Sputteringr 
process. 


U. D. C. 533.15:551.654:628.52 
J. MATHIEU 
DIFFUSION IN STRATIED BOUNDARY LAYERS 
Técnica N.º 424 — L -— 3-1975, p. 221 a 236 


The structure of the atmospheric boundary layer is studied 
2 means of a edeviation flow» from the standard atmos- 
phere, 

Mean velocity and temperature profiles are stablished for 
several stability conditions. 

Turbulent diffusion is studied both through a gradient trans- 
port approach and through a statistical theory. Difficulties 
due to the non linearity of the phenomena are analysed. It is 
shown how a destabilizing thermal action can speed up 
turbulent diffusion. 

Measurements in wind tunnel and lts limitations and advan- 
tages are compared with «in situs measurements. Difficulties 
in transpósig the wind tunnel measurements to full scale are 
mainly due to the small wind tunmel time and length scales 
and to the difficulty of simulating fairly unstable flows. 
Several ways of improving the measurements are presented. 
Economy and convenience of wind tumnel measurements 
justify its utilization and development even if our purpose is 
to seek tendencies or global effects. 


U. D. C. 621.3:517.512.2 


LUIS HENRIQUE MARTINS BORGES DE ALMEIDA 
o 
GUILHERME DINIZ MORENO DA SILVA ARROZ 


AN EXPERIMENTAL STUDY OF THE CEPSTRUM 
APPLICATION TO ECHOS AND PERIODICITIES 
DE TECTION 


Técnica N.º 424 — L— 3-1975, p. 167 a 184 


This paper, wich is a trainee report of a work done in the 
Centro de Estudos de Electrónica, shortly describes, with 
examples, the Cepstrum technique outlining its usefuliness 
in the deteccion of echoes and periodicities, and in defining 
the exact epoch they occur in. Apart from a brief theoretical 
analysis of the technique some simulations in a computer 
and some real examples are treated. 


U. D. €. 624.072.23 
J. S. BRAZÃO FARINHA 
CONTINUOUS PARTIAL FIXED END BEAMS 
Técnica N.º 424 — L-— 31975, p. 185 a 212 
The.general analysis of continuous simmetrical and partialy 


fixed end beams is developed by Cross method. Some tables 
of trapezoidal charges and rectangular charges, are presented. 
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Sociedade Geológica de Portugal 
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Instituto dos Produtos Florestais 
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Czechoslovak Academy of Sciences 
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— CENIM 
Centro Nacional de Investigaciones Metalurgicas 


— (JOURNAL OF) CHEMICAL EDUCATION 
Division of Chemical Education of the American 
Chemical Society 


— COURRIER CERN 


— CIMENTO E CONCRETO 
Associação Brasileira de Cimento Portland 


-— COMUNICACIONES ELÉCTRICAS 
International Telephone and Telegraph Corporation 


— DYNA 
Asociacion de Engenieros Industriales da Espana 


— ECONOMIC REVIEW 
Câmara de Comércio Luso-Britânica 


— EL ESTANO Y SUS APLICACIONES 
Tin Research Institute 


— ELECTRICITE DE FRANCE 
Direction des Etudes et Recherches 


— ELECTRICIDADE 
Empresa Editorial Electrotécnica Edel, Lda. 


— FUNDIÇÃO 
Associação Portuguesa de Fundição 


— FUNDEXPORT 
Fundo de Fomento de Exportação 


— GEOTÉCNICA 
Revista da Sociedade Portuguesa de Geotécnica 


— INDÚSTRIA DE MOÇAMBIQUE 
Associação Industrial de Moçambique 


— ON 
Sindicato Nacional de Indústrias Químicas 


— INDÚSTRIA PORTUGUESA 
Associação Industrial Portuguesa 


— LA HOUILLE BLANCHE 


— METALLURGIA AND METAL FORMING 
National Association of Drop Forgers and Stampers 


— MACHINE TOOL REVIEW 


— MARCONI! INSTRUMENTATION 
Marconi Instruments Limited 
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ONU (CHRONIQUE MENSUELLE) 
Service de L'Information de L'Organisation des 
Nations Unies 


O PROPULSOR 
Sindicato Nacional dos Oficiais Maquinistas da Mari- 
nha Mercante 


- PROCEEDINGS 


Mowinmlijme Nederlandse Amademie Van Wetens- 
chappen 


PLÁSTICOS 

Grémio Nacional dos Industriais de Composição e 
Transformação de Matérias Plásticas 

RAIL INTERNATIONAL 


Association Internationale du Congrés des Chemins 
da fer et de L'Union International des Chemins de fer 


REVUE D'IENI 


REVISTA DE OBRAS PUBLICAS 


- REVUE LANDIS E GYR 


Candis e Gyr sa, Service Publicité et Promotion 


REM 
Revista da Escola de Minas 


- LA RECHERCHE AEROSPATIALE 


Office National D'Etudes et de Recherches Aéros- 
patiale 


- REVISTA DE METALURGIA 


CENIN — Centro Industrial de Investigação Metalúr- 
gicas 


- REVISTA DA FACULDADE DE CIÊNCIAS DE LISBOA 


Universidade de Lisboa 
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PUBLICAÇÕES NÃO PERIÓDICAS 


- PLEUROLICINE RODENTS (GEOMYOIDEA) OF THE 


JOHN DAY FORMATION, OREGON AND THEIR 
RELATIONSHIPS TO TAXA FROM THE EARLY AND 
MIDDLE MIOCENE, SOUTH DAMOTA 

University of California 


- EARLY PLIOCENE MARINE CLIMATE AND ENVI- 


RONMENT OF THE EASTERN VENTURA BASIN, 
SOUTHERN CALIFORNIA 
University of California 


-A DELPHINID MAMMALIA, CETACEA FROM THE 


MIOCENE OF PALOS VERDES HILLS, CALIFORNIA 
University of California 


- CHIROPTERA AND DERMOPTERA OF THE FRENCH 


EARLY EOCENE 
University of California 


— OLIGOCENE BIOSTRATIGRAPHY SAN JOAQUIN 


VALLEY, CALIFORNIA 
University of California 


— MEANS OF CONTROLLING LITTORAL DRIFT TO 


PROTECT BEACHES, DUNES, ESTUARIES AND 
HARBOUR ENTRANCES 

ESTABLISHMENT OF ARTIFICIAL BEACHES 
Laboratório Nacional de Engenharia Civil 


— À HIDRÁULICA DOS ESCOAMENTOS COM FUNDO 


MÓVEL, NATUREZA, EVOLUÇÃO E PROSPECTIVA 
DE ALGUNS PROBLEMAS 
Laboratório Nacional de Engenharia Civil 


— SITUAÇÃO ECONÓMICA E FINANCEIRA DE AN- 


GOLA 
Banco de Angola — Relatório e Contas 1973 


- CONTRIBUIÇÃO PARA O ESTUDO ESTRUTURAL 


DE UMA CENTRAL HIDROELEÉECTRICA ASSENTE 
SOBRE FUNDAÇÃO DEFORMÁVEL 
Laboratório Nacional de Engenharia Civil 


— BEHAVIOUR OF A LARGE DAM BUILT ON A VERY 


DEFORMABLE FOUNDATION 
Laboratório Nacional de Engenharia Civil 


— THE INFLUENCE OF WATER TEMPERATURE ON 


THE ROUGHNESS OF ALLUVIAL FLOWS 
Laboratório Nacional da Engenharia Civil 


— ALUMÍNIO ANODIZADO PARA A CONSTRUÇÃO 


CIVIL 
Laboratório Nacional de Engenharia Civil 


— ETUDE EXPERIMENTALE DES EJECTEURS CONVER- 


GENTS TRONCONIQUES 
Office National D'Estudes et Recherches Aéros- 
patiales 


- PROPAGATION D'ONDES ACOUSTIQUES DANS UN 


CONDUIT EN PRESENCE D'ECOULEMENT 
Óffice National D'Estudes et Recherches Aéros- 
patiales 


— UN NOUVEAU DISPOTITIF D'EXCITATION VIBRA- 


TOIRE DE STRUCTURES MÉCANIQUES 
Office National D'Estudes et de Recherches Aéros- 
patiales 


- INFLUENCE DES CONDITIONS D'ENVIRONNMENT 


SUR UN RESONATEUR A QUARTZ 
Office National D'Estudes et Recherches Aéros- 
patiales 


- SYNCHRONISATION D'OSCILLATEURS ULTRA-STA- 


BLES AU MOYEN D'UN FAISCEAU LASER 
Office National D'Estudes et de Recherches Aéros- 
patiales 
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